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We consider here the problem of a "central spin", with spin quantum number S S> 1, interacting 
with a set of microscopic spins. Interactions between the microscopic spins are ignored. This model 
describes magnetic grains or magnetic macromolecules (ferromagnetically or antiferromagnetically 
ordered) interacting with nuclear spins and with any surrounding paramagnetic electronic spins. It 
has also been used to describe Si : P near the metal-insulator transition, and quantum spin glasses. 

We investigate this model in zero external field. We first set up the general formalism required to 
analyse problems in which a macroscopic degree of freedom (here, a central spin) interacts with a 
set of microscopic spins. This is done by reducing the model to an effective low-energy Hamiltonian, 
and then calculating the correlation function {S{t)S{0)) using instanton methods. Three physical 
effects come into play; we call these "topological decoherence" (coming from the phase random- 
ization by the spin environment), "orthogonality blocking" (coming from mismatch between initial 
and final environmental states) and "degeneracy blocking" (whereby the spin environment destroys 
degeneracy between the initial and final central spin states). In this paper we consider the possible 
coherent motion of the central spin. We find that all 3 of the above mechanisms suppress coherence. 
Orthogonality blocking and topological decoherence destroy phase coherence, and degeneracy block- 
ing prevents most central spins from tunneling at all. The general solution for {S{t)S{0)) is given 
for the unbiased case, for all relevant values of the couplings. Under certain conditions, nuclear spin 
diffusion also plays a role. 

The results are then used to calculate the spectral absorption x"i^) for 2 experimental systems, 
viz., TbFes grains, and ferritin molecules. For TbFez (as for most systems), coherence is completely 
destroyed. However ferritin is very unusual, in that only the degeneracy blocking mechanism is 
effective - consequently a resonance may still exist, but with a peculiar lineshape. 

Our results have general implications for the observation of mesoscopic and macroscopic quantum 
coherence, and for the foundations of quantum mechanics. They show that a spin environment 
usually has a very destructive effect on coherence, which cannot be understood using the conventional 
"oscillator bath" models of quantum environments. 



I. INTRODUCTION 



In this paper we solve a model problem which, we believe, should be of interest to both condensed matter physicists 
working on quantum spin problems, and to physicists outside this field who are interested in the foundations of 
quantum mechanics (particularly the " measurement problem" , and the physics of measurements and " decoherence" ) . 
It should also be of interest to physicists and chemists working on quantum dissipation and models of the environment, 
for one of the central points in what follows is the elucidation of the properties of a new kind of quantum environment, 
quite different from the usual "oscillator bath" models. 

The model we discuss was originally introduced to deal with the effects of environmental spins on the coherent 
tunneling of magnetic grain magnetisation fll-^ . It consists of a central " giant spin" , having spin quantum number 
S' ^ 1, coupled to a set {cjk} of "microscopic spins"; the microscopic spins are not coupled to each other (Fig.l). The 
model is then completely described by a Hamiltonian Ho{S) for the central spin, a Hamiltonian Henv{{'^k\) for the 
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microscopic spins, and a coupling Hint{S, {f?fc}) between the two. One can also introduce an external field Ho which 
couples to these spins - the effect of external fields will be considered in detail in a 2nd paper. 

In most practical applications of this "central spin model", the environmental spins will be either nuclear spins 
(which may be either inside the object carrying the giant spin, or outside it, in some substrate, or solvent, or 
surrounding matrix) or else paramagnetic electronic " spin impurities" , which may also be inside or outside the giant 
spin. The central spin may be a magnetic grain Q or a magnetic macromolecule such as ferritin j|] or, on a smaller 
scale Mni20i2\ or it may be one of the large superparamagnetic "spin clusters" which are believed to exist in many 
disordered magnets at low temperature, such as Si : P near the metal-insulator transition ||^, or "giant magnetic 
polarons" Similar spin clusters exist in "quantum spin glasses" 

The crucial point which has come out of our study of this model, and which we believe is of great importance for 
both quantum magnets and for discussions of decoherence and the theory of measurement, is as follows. The effect 
of microscopic spins on the quantum coherence of the central spin is akin to the effect of an invasion of microscopic 
viruses on a host bacterium, i.e., in most cases complete destruction! Environmental spins act as a far more potent 
suppressor of quantum coherence than any other kind of quantum environment hitherto discussed - in particular, the 
conventional bosonic oscillator baths ||l(]|-|l2|] seem relatively benign. 

Now although the spin environment we consider constitutes an example of an " unconventional environment" ref . , 
a moment's thought shows that such spin environments will be the rule rather than the exception, when considering 
macroscopic or mesoscopic quantum coherence. This is because most large quantum objects will have at least a 
coupling to their own nuclei. Thus in the oft-cited example of the superconducting SQUID ring fl^ , there is an 
electromagnetic coupling of the tunneling flux coordinate (and the associated ring supercurrent) to all of the nuclear 
spins within a penetration depth of the surface (either in the junction or in the ring), as well as to any paramagnetic 
impurity spins that happen to be present [ p^ . While the effect of these interactions on tunneling may not be serious, 
except under special circumstances their effect on coherence (MQC) will be quite drastic. 

This means that our results are not only of practical importance for any experiments which attempt to test quantum 
coherence on the mesoscopic or macroscopic scale. They also show why coherence is so easily destroyed at the 
mesoscopic or macroscopic level (in particular, in measuring devices). This point is not obvious in the previous 
discussions of, e.g., MQC in the conventional models, in which bosonic oscillators models of the environment are used 
- in fact, it is predicted, using such models, that MQC ought to be visible in SQUID rings p^ , p^ and in magnetic 
grains Q,^, using presently achievable values of the classical dissipation in these systems. Whilst no observation 
of MQC has yet been claimed in SQUIDs, it has been claimed that MQC has been detected in samples of ferritin 
magnetic macromolecules [^0. Our central spin model can be directly applied to this case - as we shall see here, 
the case of ferritin is very unusual, in that most of the decoherence mechanisms we shall discuss are inoperative. In 
fact only one mechanism, which we call " degeneracy blocking" , operates on ferritin, leading to an unusual lineshape 
which is not obviously inconsistent with the experiments (the ferritin experiments have also been criticized on other 
grounds, mostly related to the power absorption and the behaviour as a function of field |lj]). In section^ we will 
give a quantitative discussion of spin environment effects on ferritin, and also TbFe^ grains. 

Our general thesis is that the almost universal presence of spin environments is going to make it far harder to see 
quantum coherence on the macroscopic scale than has previously been realized. Even on the mesoscopic scale it will be 
difficult. There are general strategies that one may adopt in looking for suitable candidates for MQC (or mesoscopic 
coherence). For example, it helps to choose a quantum soliton as one's macroscopic quantum object [ p"9|j2^ , since 
quantum solitons, by definition, do not couple linearly to other excitations of their own field. An example of current 
interest concerns the quantum tunneling of macroscopic magnetic domain walls - recent experiments have apparently 
discovered this in Ni wires |2l| , and it will be interesting to see how experiment compares with theory in this case. 
It is also useful to make the quantum object neutral - this removes all electromagnetic couplings to the environment 
(although these are often not important, since they are infra-red weak [^). However the crucial step is clearly one 
of getting rid of any environmental spins - this suggests at the very least a system for which the nuclear spin is zero. 
There is in fact an almost ideal candidate [l4| , viz., superfluid "^He. Here the macroscopic quantum object (whose size 
is defined by the geometry of the container) is a vortex (a special kind of topological quantum soliton) . The system is 
neutral, and the nuclear spin is zero; moreover, the system can be purified to a quite extraordinary degree p^ , and 
the experimental cell walls coated with solid ^He to remove any influence from impurities, bumps, etc. In fact the 
only couplings left are to phonons; in common with other quantum vortex problems, such as magnetic vortices [ pO[ , 
this coupling has both a local and a non-local component. However, just as with photons, that part of this coupling 
which is responsible for decoherence is infra-red weak (both of these couplings can be described by a Caldeira-Leggett 
spectral function J(w) ~ uj'^). 

In this paper we will not be able to discuss all of these issues, concerning the physics of spin environments. In 
particular, we intend to discuss elsewhere some of the consequences for both mesoscopic quantum devices and quantum 
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coherence experiments. These apphcations of the theory are at the forefront of present attempts by physicists to push 
the "F.A.P.P. barrier", between the quantum world and the classical world |2^, ever further into the mesoscopic and 
even macroscopic domain. The implications for the measurement problem will also be dealt with elsewhere. 

In what follows we begin (in section ||) by describing the central spin model, and giving some discussion of the 
relevant physical coupling energies in it, with the aid of physical examples (ranging from magnetic grains to "spin 
clusters" in Si:P). In section |ll| we show how an initial microscopic Hamiltonian, describing the giant central spin 
and its microscopic environmental satellites, can be reduced to a low-energy effective Hamiltonian description 
The particular way this reduction is done, and the terms that are left at the end, is partly a matter of theoretical 
choice. We have done it in a way which displays most conveniently the important physical processes that operate at 
low energies. To help the reader who finds this reduction a little abstract, we sketch how this reduction is done in 
practice for a real system. 

In section IV we proceed to the discussion of 3 solvable limits of our model Hamiltonian |^,||. Each of these 
limits brings out a different physical effect. The simplest of these to understand is "topological decoherence" , in 
which the topological phase of the spin environment becomes entangled with the coherent phase of the giant spin, and 
destroys the phase coherence of the latter. This decoherence mechanism is completely novel, and has many amusing 
consequencies - it leads to decoherence without dissipation or energy exchange of any kind, and the decay of the 
correlation functions of the central spin, caused by it, is most unusual ^,||. The second effect |^ is of "orthogonality 
blocking"; it arises when there is a mis-match between the final-state spin wave- functions in the environment, and 
their initial states. Although this effect is reminiscent of the "orthogonality catastrophe" first discussed by Anderson 
, which is at the heart of decoherence in the conventional oscillator bath models , it differs in important details 
(for example, instead of exponential suppression of tunneling rates, one gets power low suppression; and the effect 
on the correlation functions of the central spin S is most unusual). Finally we study the limit appropriate to pure 
"degeneracy blocking" [§|J|], which arises because the coupling between the central spin and the environmental spins 
lifts that degeneracy, between initial and final states of the central spin, which is essential for coherence. 

In this paper we will largely concentrate on the coherence properties of the central spin, rather then its relaxation 
from any particular state. As we shall see, this usually requires consideration only of processes in which the total 
polarization of the environmental spin bath is unchanged when 5^ fiips. Other processes necessarily destroy coherence. 

Having thoroughly understood the 3 limits, we proceed in section ^ to the discussion of the "generic case", i.e., 
having arbitrary values of the parameters in the effective Hamiltonian, so that all 3 mechanisms are operating. At 
this point we also include the weak nuclear spin diffusion effects, which can also play a role in decoherence. 

In section VI, we return to two of the physical examples described in section |l|.C (i.e., TbFe^ and ferritin), and 
show how their behaviour is determined by the solutions derived in sections IV and ^ In the case of ThFe^., which 
is fairly typical, we find a catastrophic destruction of coherence - no trace of any resonance is left. On the other hand 
ferritin has very few nuclear spins, and one finds that only degeneracy blocking play any role in the final behaviour. 
As a result there is a significant reduction in the spectral weight of the resonance, but the lineshape still shows a 
strong but highly asymmetric peak. Finally in section VII, we summarize our results. 

An important omission from this paper is the effect of an external field. This field will strongly affect the giant spin 
dynamics (with usually a much lesser affect on the microscopic spins). Moreover, dissipation can play an essential role 
in the problem, whereas in the unbiased case we consider, decoherence arises without dissipation. A full discussion of 
the relaxation of {S{t)S{Qi)) (as opposed to coherence) also demands that we consider incoherent transitions involving 
a change in the spin bath polarization - moreover, since the relaxation of both the central spin and the spin bath is 
determined by coupling to phonons and possibly electrons, these must also be included. In view of the experimental 
importance of the field effects, we have reserved discussion of them to a second paper. 

For practical applications of our results, it is important to note that there will also be an influence on coherence 
coming from phonons or electrons. These can interfere in interesting ways with the effect of the environmental spins 
|p5| . We hope to address this aspect in more detail at a later time. 

Many of the results appearing in sections |] , III and IV of the present paper have appeared previously in short 
communications [p]-p|j2^]. Here we give detailed derivations for the first time. The results in section^ (the generic 
case) are new here, as is the detailed analysis for the various applications. Some of the mathematical derivations are 
quite lengthy, and are relegated to Appendices. 



II. THE MODEL, AND SOME SIMPLE RESULTS 

In this section we give a physical discussion of the central spin problem, with reference to a number of different 
examples, including magnetic grains, "spin clusters", etc. The various physical interactions in these problems are 
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described and quantified. 



A. The central spin 

A typical central spin is made up of a very large number of microscopic spins, whose motion is locked together 
by nearest-neighbour exchange couplings (or, in the case of superparamagnetic spin complexes, by indirect spin-spin 
couplings such as the RKKY interaction). The exchange coupling can be very large, as much as 1 — 2eV; this is much 
larger than other individual spin energy scales such as those coming from anisotropy, so that in principle one can have 
large monodomain magnetic particles, with many spins (up to perhaps 10*, depending on the particular systems) 
lined up to form a " giant spin" . Another more subtle example is the antiferromagnetic grain, in which the sign of 
the nearest- neighbour exchange is positive, so that nearest neighbours are anti-parallel, and one obtains a " giant Neel 
vector" l^?!] . The superparamagnetic spin complexes are more delicate still, and typically arise at low temperatures 
when paramagnetic electronic spins in a disordered magnet "lock together" via indirect exchange, RKKY interactions, 
etc. 

Now it is clear that a realistic Hamiltonian for such a giant spin will be very complex. Let us consider first the 
simplest case of a ferromagnetic grain. An enormous literature exists on the properties of such grains, ranging from 
simple elemental grains, such as Ni, Fe, Tb, or Ho, to more complex magnetic compounds such as FeCo^, Fe20^ 
(7-haematite), Mni20i2 and its derivatives, right up to organic magnetic molecules. With present techniques it is 
possible to prepare many grains of very uniform size ||2^ ] (or else rely on nature for the preparation of magnetic 
macromolecules) . 

An obvious first step in analysing such a system is to treat the entire spin complex as a rigid quantum rotator S, 
with dynamics governed by the single-spin anisotropy field in the particle/grain; i.e., we write a central spin "bare 
Hamiltonian" Ho{S), with | 5 |= 5 a constant. Early studies of the quantum dynamics of this model were stymied 
by the constraint of constant | 5* | , which implies a kinetic term in the Lagrangian which is linear in time derivatives 
- this leads to well-known difficulties when one tries to apply path- integral methods to the problem |2S[| . A number 
of crucial advances allowed the solution of the problem, most particularly the proper treatment of the kinetic term, 
which leads to the " Haldane phase" , and also a correct formulation of WKB theory for the tunneling dynamics 
of S, which yields the correct semi-classical solution for tunneling of an angular momentum through the barrier . 
This semi-classical solution was also later derived using the instanton calculus (although there still seem to be 

differences between the 2 methods, concerning the tunneling prefactor |^l| ). 

As a simple example of the bare Hamiltonian Ho{S) we may take (we write S = Ss, where S* =| S |) 

HoiS) = -iq.sl+K^sl, (2.1) 

which has two degenerate classical minima at Sz — ±S; the semi-classical paths between these 2 minima move on or 
near the easy XZ-plane (note that iiTy and K± are proportional to S). Notice that if K± = 0, there is no tunneling 
because Sz is conserved; this point is obviously connected with the fact that we are dealing with angular momenta and 
not particles, and caused som e confusion in the literature before 1986. Application of a magnetic fi eld adds a term 



—M ■ Ho = —^S ■ Ho to (p. If ); if Ho is directed along z, this biases the symmetric 2-well problem in (2A), but if it is 
applied along x, it lowers the barrier height and displaces the degenerate minima towards each other in the XZ-plane. 
Application of Ho along y distorts the semi-classical path and lowers the barrier. A further effect of an applied Ho 
is mostly nicely seen in the instanton language; since the 2 possible paths between the degenerate minima involve 
opposite Haldane topological phase |^3| , this phase can be changed by an external field p4[] , causing the tunneling 
splitting to change and even oscillate. 

In fact the original papers of van Hemmen and Suto analyzed a very general class of central spin bare Hamiltonian, 
given by 

n 

Ho{S) = -i^ii I Sz I' -1/2^ ^Kr (s; + sL) , (2.2) 

r=l 

and a considerable number of more recent papers have concerned themselves with re-analyzing special cases of ( [j.2| ). 

I). One obvious omission is the infra-red weak coupling of 



Consider now what has been left out of (2.1) and (2.2 



S to photons (via the magnetic dipolar interaction ||19,^|) and to phonons. Another omission is coupling to higher 
excited states of the particle, involving spin fiips of spins inside the grain (internal magnons); the neglect of these is 
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usually justified by appeal to the large value | J | of spin exchange y,^, in models like equation (|2.lD. Incorporation 
of such processes implies relaxation of the constraint | S |=constant. Another omission is that of surface excitations 
of the grain - in any real grain, there can be low-energy magnetic excitations involving "loose spins", i.e., spins which, 
because of imperfect preparation or the inevitable defects, dislocations, etc., have a coupling energy <C J to their 
neighbours. The effect of these is hard to quantify, but will be included below in our description of environmental 
spin effects. 

Some of the above additional physics has been discussed already in the literature. The effect of phonons has been 
analyzed in the usual Caldeira-Leggett instanton scheme by Garg & Kim and Chudnovsky . These calculations 
indicate very small corrections to tunneling rates coming from phonons, which is not surprising given the infra-red 
weakness of phonon environments. Notice however that these calculations ignore barrier fluctuation effects [^7[-p9[. 
From now on we will ignore phonon effects since they are so small. For the same reason we will ignore the even 
smaller dissipative effects of photons (calculated in ref. [^). However we note again that in quantitative discussion 
of relaxation experiments, as opposed to coherence experiments it will be necessary to include phonon effects. 

We shall also ignore coupling of S to electrons. Actually, in contradistinction to domain walls |4^,^l[], electrons 
have a very powerful effect on the dynamics of S. Here we simply assume the grain is insulating. 



We now turn to the really important omission from (2.1) and (p.2|), viz, the coupling of S to other spin degrees of 



freedom - as already noted, these will include nuclear spins and paramagnetic electronic impurities, both inside and 
outside the grain. It will not include the low-energy excitations of coupled spins inside or outside the grain - these are 
generally describable in terms of a bath of oscillators (e.g. magnons) with a very weak super-Ohmic coupling to S. 
Thus it is crucial that the environmental spins be weakly coupled to each other - then the variety of effects we study 
in this paper become possible. 

The most ubiquitous of such environments is the nuclear spin environment - this is well understood. A nucleus with 
a finite spin /, and spin moment — 7jv?i/, will interact both locally with the currents of electron clouds at the same 
ionic site, via the contact hyperfine interaction, and also non-locally with other ions via the dipolar field. There may 
also be other residual interaction such as "transferred" hyperfine interactions (as in, e.g., Fe^^ in MnF2 or C0F2), 
or quadrupolar couplings if / > 1/2, or nuclear spin-phonon interaction. Finally, one may have interactions between 
the different nuclear spins, via, e.g., spin waves (the Suhl-Nakamura interaction), or dipolar couplings - these are 
very weak (Suhl-Nakamura interactions are < 10"^ K, and nuclear dipole-dipole interactions ^ 10^^ K) but, it turns 
out, can play a role in the dynamics of S; this is discussed in Section 0E. The most important interactions are the 
contact hyperfine interaction and the dipolar interaction between the nuclei and electrons. Thus the nuclear moments 
in nonmagnetic hosts (e.g., in H2O, or in CuClz ■ H2O) will be described by a Hamiltonian 

Hi = -iNhl- |go-7e?^E;|((^^(^)) -3 '^"^'^ '|^^^^^^ )| - (2-3) 

where Ho is an applied field, ?i7e = /isffe, where is the electronic g-factor, sums over all ionic sites and {Sj{T)) 
is the temperature-dependent expectation value of spin polarization at the sites (for weak fields, {Sj(T)) = XjHo, 
where Xj is the susceptibility tensor). On the other hand a nucleus in a magnetic ion will have a Hamiltonian 



Hr 



o ) + Ho + Hdip^ , (2.4) 



//.,, = -7e?.E3(4-3^1^). (2.5) 

The strength of these interactions is measured directly in NMR experiments; the hyperfine coupling range from 
1 — iMHz (i.e., 5 — 15 X \Q~^K) for protons in _ff^, up to values greater than 5000MiJz (0.25 X) for some rare-earth 
magnetic nuclei (e.g., Tb^^^, Dy^^^); these latter correspond to local fields acting on the nuclei as high as 500 TesZa, 
and come overwhelmingly from the contact interaction. Thus when a central spin rotates, with all internal electronic 
spins rotating in unison, the contact interaction tries to force the nuclei on the magnetic sites to follow. If the hyperfine 
interaction strength is ujo, then it is clear that the ratio Wo/f^o, with the timescale of central spin rotation, is 
going to be crucial to the nuclear spin dynamics. If the electronic spin frequencies flo 3> Wo, then the nuclear spins will 
experience a "sudden" perturbation, and few of them will follow the central spin; conversely, if tOo 3> the nuclear 
spins will follow adiabatically. In fact, as we shall see in more detail later in the paper, the ratio oJo/^o is usually 
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somewhere between 0.01 — 0.05 (although for rare earth nuclei like Tb^^^ or Dy^^^ it can be ^ 0{1), depending on 
the host). On the other hand for nuclei in non-magnetic ions, tUo/^o may be much smaller, and will of course depend 
strongly on the host (in a magnetic host, dipolar interactions from permanent moments, and also possible transfer 
hyperfine interactions —"/nHI ■ Aoj ■ Sj from nearby magnetic moments Sj, can greatly increase cOo', for example, lUo 
for in MnFa is 160 MHz). 

A large central spin will also interact significantly with nuclear spins in a surrounding medium, such as a substrate 
or solvent, via the dipolar interaction generated by the central spin dipolar field, i.e., 

= ~n^le,N E ( E ;i - ,fAjA^ ) , (2.6) 

where is the distance between some nuclear moments It in the surrounding matrix and a local moment Sj which 
is incorporated in the central spin. This interaction is not negligible; for a central spin with S = 10^, the dipolar 
coupling to a nucleus at a distance of 10"^A is already ~ 1 MHz, rising to 30 MHz at a distance of 300A; and there 
is clearly a very large number of nuclear spins within lOOOA of the central spin! 

This brings us to another potentially very serious source of decoherence for the central spin, viz., paramagnetic 
impurity spins in the surrounding matrix. The Hamiltonian has the same form as above except that 7jv^i is replaced 
by 765*;, describing an impurity electronic spin Si at some distance rn from the central spin. In a very pure surrounding 
matrix this is a problem because 7^ ~ 2 x 10'^7Ar; for the central spin with S = 10^, this dipolar coupling is still 
> 2 MHz for a paramagnetic impurity 10* A distant! A sphere of this radius will contain over 10^^ ions, so unless we 
have a matrix of, e.g., solid ultrapure ^He (or perhaps very high grade Si), there will be many such impurities with 
appreciable coupling to S. 

Finally, as we approach the surface of the giant spin, we may expect a number of electronic spins to have "weak 
exchange" coupling to the giant spin. This may arise from imperfections or defects or even dirt at the surface (or 
near it, as in the coating around ferritin molecules). Typically we expect some kind of " superexchange" coupling to 
apply, with energy anywhere between 10 — 1000 K (i.e., roughly in the range 10^" — lO^'^ Hz). We shall refer to such 
spins in the spin environment as " loose spins" . How many of them there are will depend very much on the materials 
and material preparation used. 

A somewhat similar picture may well apply to Si : P near the metal-insulator transition, where it is believed that 
the ESR lineshape ji^] may well be explained in terms of a Hamiltonian 

N 

H,ff^AS-J2ckIk+Ho-S, (2.7) 

k=l 

in which a central spin "superparamagnetic cluster" interacts with nuclear spins. This may well be so, provided there 
are insignificant decohering interactions between the spin complexes. We hope at some time to give a detailed analysis 
of this case, and also of similar quantum spin glass systems. 

We may summarize this survey of the various components of the spin environment surrounding a ferromagnetic 
grain by the diagram in Fig. 2. It is important to notice how completely the range of frequencies from 10^ Hz 
downwards is covered by the various coupling strengths; and even above 10^ Hz, for very large grains and/or including 
"loose spins". 

B. Some examples 

The reader of this paper who is not a specialist in magnetism may also find it helpful to get a feel for specific 
numbers, for different materials in which different interactions are important. Thus in this sub-section we will give 
some "vital statistics" for a number of real magnetic systems. 

Before doing so, it will be useful to glance at Table I, which lists some nuclear hyperfine coupling energies, in 
various host systems. For those who like to think of energy in temperature units, note that 1 K ^ 21.6 GHz. It will 
be noticed that the couplings range over a factor ^ 10'^, and that even the nuclei from non-magnetic ions can have 
reasonably large couplings, via transferred hyperfine interaction {F^^ in the Antiferromagnet (AFM) MnF2 being a 
good example). For the transition metal Ferromagnets (FM), oj^ ranges from 50 — 500 MHz, but the rare-earth FM's 
can have Uj- an order of magnitude greater. Where more then one value of is given, there is either a quadrupolar 
splitting of the fines (e.g. Tb, Dy) or nuclei are at inequivalent lattice sites (as in YIG). 
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The internuclear dipole interactions are way down from this (typically ~ 0{lkHz)); and even the indirect Shul- 
Nakamura interactions have energy scale ~ 100 -200kHz (see Table ||). Of course these interactions are always there, 
and they lead to an inevitable minimum spreading Su)k in the couplings; but Scuk/uik is clearly very small (10~^ — 10""^) 
and so our assumption of independent nuclei is very accurate (note that the temperature T will always be far larger 
than this linewidth Sujk)- However we notice that Scok is not necessarily small compared to A,,; this will be important 
for the magnetic relaxation at long times, but is irrelevant for coherence, as we shall see later. 

We have already seen (cf. Fig. 2) how the coupling to external nuclei and external electronic spins is also important. 
As a rough guide, this coupling can be estimated as ujk ^ 10^ S/r^ Hz, where S is the spin quantum number of 
the central spin, and r is the distance of a nuclear spin from S in A. For an electronic spin, 2 x 10^^ S/r^ Hz 

(obviously the prefactor is modified depending on the spin quantum numbers and ^-factors of the environmental spin). 

Consider now a few experimental systems, by way of example: 

(a) MnF2 grain: This is an AFM; however it will have an uncompensated excess spin moment which may couple 
to external nuclei or electronic spins, depending on the host. The size of fio will depend on many factors for a grain, 
mostly uncontrollable strain fields - in principle one expects fio to range from ~ 0.1 to IK for such a system. 
There are two hyperfine couplings; Wfc = IQOMHz {F^^) and 680 MiJz (Mn^^). If fio — 0.1 K this gives values of 
ctk = {t^ /2)u)k/^o equal to 0.125 and ~ 0.5, respectively; if fio = 1 if, they are ten times smaller. 

(b) YIG grain: One can make very pure grains of Yttrium Iron Garnet. Since it is a soft magnet, a well-prepared 
system would have Vto < 0.1 K, even perhaps as low as 0.01 K . In the former case, the Fe^"^ hyperfine coupling has 
ak ~ 0.05; in the latter case ak ^ 0.5. Note that only 2.19% of Fe nuclei have a moment. 

(c) Ferritin : Ferritin is a magnetic macromolecule made in all eukaryotic cells; a typical molecule contains 4500 Fe 
ions in a ferritin structure, which are antiferromagnetically ordered, although there is an excess uncompensated 
moment of size somewhere between 200 and QOOfXB- The 2.19% of Fe^^ ions (roughly 100 of them) have Uk = GAMHz. 



In the experiments of Awschalom et al. 17|, an anisotropy field ^ 1.72 Tesla was found, indicating f2o ~ 2 i^T, which 
indicates a value for ak ~ 2.2 x 10~^. 

(d) TbFes grain: Here JIq is typically 3 K; this is very high indeed. However so are the ojk since Tb is a rare-earth 
magnet. With the quadrupolar split hyperfine lines, one finds ak ~ 0.08, 0.07, 0.05 for the Tb nuclei. The 2% of Fe^"^ 
nuclei will have ak some 50 times smaller than this; the result of this for coherence of S will be, in effect, to simulate 
a spread Sujk in the Wfc, rather as though we had a spread Sujk ^ 50 MHz in the system. 

From these examples we see that typical values of ak are between 10^"^ and 0(1), usually something like 0.05. 
Moreover the presence of other nuclear species such as H^ (in H2O, as in ferritin) will cause a small spread in the 
ak- Finally, of course, there will be a much larger spread in the ujk coming from spins outside the grain or molecule. 

In section VI we shall return to examine two of these examples, namely ThFe^ and ferritin, in greater detail. 



III. THE EFFECTIVE HAMILTONIAN 



In this section we make the crucial step of going from "bare Hamiltonian" like those in Section to a low energy 
" effective Hamiltonian" . This effective Hamiltonian is our starting point for calculating the dynamics of S at low T. 
It contains a variety of effective couplings, which already contain all the high-energy physics. We emphasize that in 
realistic cases one may actually calculate these renormalised couplings, starting with the original bare Hamiltonian - 
they do not have to be treated as phenomenological couplings to be determined from experiment. To illustrate this 
we actually do the truncation explicitly for a typical hyperfine interaction. 



A. Reduction to i/e// 



Let us go back to our "bare" Hamiltonian Ho{S) for S, before either truncation to a low energy effective Hamiltonian, 
or coupling to any environment. There are two important energy scales in Ho{S), viz, f2o and Aq; the former describes 
the "fast" motion of the central spin in the effective potential in Ho{S) (small oscillations near the potential minima, 
or the "bounce frequency" involved in tunneling through the barrier), whereas the latter characterizes the much 
slower coherence between tunneling events, since 7r/2Ao is the average time between successive tunneling events. 
These 2 time scales are seen very clearly if we show a "typical trajectory" for S{t), i.e., a typical path that would 
contribute to a path integral evaluation of the dynamics of S. This is shown in Fig. 3, for a case in which Si and S2 
are not antiparallel. There will be other even more rapid time scales which are buried in this path, corresponding to 
transitions to very high excited states of Ho{S); but their amplitude is very small. 
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The truncation of Ho{S) is standard |3^j3l|,Q. Consider, e.g., the bare Hamiltonian truncation gives 

HoiS) ^ 2A„f, cosTT^ = H,{f) , (3.1) 

where the PauU matrix flips S between the degenerate configurations Si = zS and 5*2 — —zS, the factor 2cos7rS' 
is an extression of Kramer's theorem. The bare spUtting Aq is given in terms of fio by 

Ao ~ 12oe-^° , (3.2) 

Ao^2{K\\/K^f'^S , (3.3) 

showing the exponential ratio between Aq and Q,o (more accurate expression for Aq appear in Refs. |,|^,|l|); and 
for this model 

r!o~ |(/^l|i^^)i/2 . (3.4) 



of (|3.2D 



No tice that ahhough Vto is absent from (3J_), it is nevertheless implicitly contained via the " WKB renormalization" 



Let us now consider how one should make the same truncation in the presence of a large number of environmental 
spins. In general the same argumentation will apply - we look first at how the environmental spins will affect the 
"fast" physics of the central spin (on a time scale ^), and absorb this into a new effective Hamiltonian describing 
a " renormalized instanton" for the central spin. We then deal with the "slow" physics (on time scales > il^ ^) using 
this renormalized instanton, and derive a new effective Hamiltonian for this slow, low-energy physics. 

However the same is certainly not true if we have N environmental spins. How can we handle this case? The 
method we have adopted [^,|| starts by asking how we should write the effective action to be associated with the 
instanton. In the absence of the spin environment, the transition amplitude matrix from | 5*1) and | 5*2) is given by 

= tMo exp{~Ao ± inS} - Aof,e±*"^ , (3.5) 



where Aq = floe and Ao ^ {K^^/ K±)^^'^S, in accordance with equations ( ^^ ) and (3.4); the ± refers to the 2 
different (clockwise or anticlockwise) paths between | 5*1) and | 52). Now when we add the spins, the instanton will 
act as an operator in the subspace of each environmental spin (able to cause transitions in each subspace, in line with 
our discussion in section |l|.B). The general form now must be 

(k) = Aof^e*'' exp | [P^Uk ± ^kVk] ■ ± ivb^S + 4>k + ctkfik • | , (?? = ±) . (3-6) 

in the subspace of the /c-th spin - here f,, = t± = f ± iTy, Uk, n^, and Vk are unit vectors in some direction, and 5k 
is a number. The instanton from | ^2) to | Si) is Hermitian conjugated to that in ( |3.6| ). The role of the new terms in 
( |3.6D is as follows: 

(i) The adiabatic correction 5k is due to the high-frequency (3> flo) fluctuations of ijk, which renormalize the bare 
instanton action A^. This correction is usually very small, and is of course unobservable because it does nothing but 
renormalize Aq. It increases the "moment of inertia" of S. 

(ii) The "barrier fluctuation" terms {(3kUk ± Cfc^fc) • come from fluctuations of ak which are of frequency flo 
or less. As discussed in refs. |3^ and (see pp. 53-55 of ref. [Q), the effect of these fluctuations is, amongst other 
things, to cause the potential barrier to fluctuate in height. In Section ^l| B we discuss in more detail how this leads 
to a term of this form for our problem. Here we simply note that (a) the effect of is to increase the tunneling 
rate, because barrier fluctuations occasionally lower the barrier , thereby creating a "window of opportunity" for 5' 
to tunnel; and (b) this term can be entirely absorbed into a renormalization of the "orthogonality blocking" term 
that we will have in our flnal effective Hamiltonian. In the most general case one has to allow different couplings for 
the clockwise and counter-clockwise instantons, e.g., when the environmental spin is oriented along the x direction it 
influences the two instanton trajectories in the XZ-p\a,ne differently, thus giving a non-zero value of ^k- 

(iii) The topological term akfik ■ dk + 4>k describes the dynamical effect on Sk when S flips. Its origin is two fold: 
one contribution comes from a transfer matrix (note that ak ^olk) 
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in the spin subspace of the fc-th spin, where in the strong couphng case the phase 4>k approaches the value (j)k = 
(j)^ — 7r/2; (f)^ is the adiabatic Berry phase accumulated by the nuclear spin. The other contribution {ak — ctk) is 
closely related to the barrier preparation effect, but has no analog in the standard description of the particle tunneling 
through the barrier; it arises from the influence of the environmental spins on the topological phase (see below) . 

This summarizes the mutual effects of S and {crk} during the tunneling of S. However after this tunneling is over, 
there will be a residual interaction between S and the {fffc}- This can be described as follows. Suppose the 2 relevant 
quasiclassical states of S are again oriented along 5*1 and 5*2, and that the total effective fields acting on (7^ are 
correspondingly ■f^^-' and 7[.^\ defined in units such that their interaction energy with ffk is 7fc ■ (fk- Now define the 
sum and the difference terms 

r -(1) -(2) 

(3.8) 

where the Ik and irik are unit vectors. In the truncated Hilbert space of f and the residual interaction term 

takes the form 



N N 



1 f II - - 

k=l k=l 



(^k -nik- Sk \ , (3.9) 



i.e., a term which changes when Si — > S2, and an extra term which does not. This latter term can arise in various 
ways. Most commonly it will arise because the spin-1/2 variable cjk is produced by truncating some environmental 
spins /fc, in a bare Hamiltonian, with Ik > 1/2. It will also arise if, e.g.. Si and ^2 are not exactly antiparallel (e.g., 
in an imperfect grain) , or simply because there are small stray fields acting on the ak ■ Usually to^ <C LUk (^^ least for 
nuclei in a magnetic system). 

We are now ready to write the full low-energy effective Hamiltonian for our central spin model, in the absence of 
an external field, but incorporating all the mutual effects between S and its spin environment |^,^ 

N N \\ N j_ 

Heff = 2Ao{f+ COS [$ + ^(afcnfc - i£,kVk) ■ ^k] + H.c.} + ^ ^ 4 • + ^ ^ mfe • . (3.10) 

k=l k=l k=l 

The form of this effective Hamiltonian is central to all further discussions. We have dropped the barrier fluctuation 
term (3kUk • since as already noted its effects will simply consist in a renormalization of Aq and the effective ratio 
'^k/^h whose effects appear in orthogonality blocking (see the next section). The phase $ = {nS + J2k^k), i-e., we 
absorb the Berry phase into our total value of the topological phase. Finally, Aq — Aoexp{— J^k ^k), absorbing the 
adiabatic corrections. 

Before going on, it is worth asking (a) how general is Heff, and (b) what, if anything, has been left out of it? 

Consider first the general structure of Heff. We see that in our 2^^+^-dimentional Hilbert space, it contains not 
only all possible pairwise interactions between r and each tTfc, but also a set of higher couplings, involving as many as 
N environmental spins, produced by expanding the cosine; i.e., we have couplings of form 

Hdyn ^ TiAl^p^^s...ak^a'^^d-l^al^ ■ ■ ■ . (3.11) 



Some readers might like to think of ( 3.11 ) as being an infinite set of terms i n a p erturbation expansion of the bare 
coupling to all orders, but we feel that this is mistaken - the coefficients in ( 3.1l|) are fixed precisely by the cosine. 



and in any way there is nothing perturbative about the form of Hef / - it is valid even when the coupling between the 
ak and r is very strong. 



In fact the only terms left out of Heff are direct interactions between the {ak} (i-e., ones not generated, like (3.11), 
via coupling to S). For nuclei these are the dipolar and Nakamura-Suhl interactions, and for paramagnetic impurities 
they could also be RKKY interactions between widely-separated spins. Thus the only terms left out of Heff are terms 
of the form 

Viiak}) = ^kl'^l^l' ■ (3-12) 

k^k' 

As we shall see later in this paper, these terms can also be included in our analysis. Their effect on coherence is 
usually small (although under some circumstances they can effect coherence in large grains). However in a second 
paper, dealing with relaxation of 5* in a bias, we will see that V does play a very important role. 



9 



Readers not used to the philosophy underlying the derivation of low-energy effective Hamiltonians may also suspect 
that somehow the form of H^ff depends in some way on the semiclassical analysis we use of the dynamics of S. We 
emphasize this is not so - the form of Hgf f depends in no way on any semiclassical or WKB approximation. What 
does depend on the semiclassical analysis is the exact value of each of the renormalised parameters in Hf,f f.liS were 
a microscopic spin, it would not be appropriate to use WKB theory to evaluate, e.g., Aq in Ho{t), or the {ak} in 
Heff. However when S* ^ 1, the use of WKB/instanton methods becomes extremely accurate, both in determining 



Ho{t) (cf. Rcfs. |g9|,|T-33|), and (as we shall see in the next sub-section) in determining the values of the parameters 
in Heff. 

Finally, note how differently H^ff describes the "spin bath" environment from the usual oscillator bath llo|-|l^] 
models. Not surprisingly, the physical effects of the spin bath will also turn out to be very different - from this point 
of view it is useful to think of the spin bath as an example of an " unconventional environment" , having unusual 
effects (whether it interacts with a central spin, or SQUID, or any other collective coordinate). This theme has been 
elaborated elsewhere 1l3|. 



B. Instanton Derivation of -ffe// 
As a first step we will present general arguments and the physical ideas involved in such a derivation, and proceed 



then with the more specific calculation for the Hamiltonian (2.1). 

Suppose we are interested in the transition amplitude from | S*!) to | 5*2) in the WKB approximation for the central 
spin only, ignoring for the moment the initial and final states of the environment. Without any interaction between 
the S and {Bk} we have, of course, 

= Ta^o exp{-A„ ± i^S] , (a = ±) , (3.13) 
Now if we calculate the transition amplitude with all the interactions included, according to 

I )}) = {S, I e-^!dr'[H^(r')+H^^Ar')\ | 5?^^ | |^(»)|^ ^ (3^14) 

it will obviously depend on the initial state of the {(t^}, and, when projected on the particular final state of the 
environmental spins, it will depend on the spin flips in the environment. Thus Mm in the most general way we have 



to consider the instanton from | Si) to | S2) as an operator in the subspace of {Sk}- Eq.(|3.6D is the most general 
form for this operator one can write down no matter what are the interactions in the system (the only omission is 
the simplified form for the adiabatic contribution 5^ which in a more general form should read (5^ ± (5^, applying to 
the case when the central spin is in an external magnetic field so that the adiabatic contribution is asymmetric for 
the clockwise and anticlockwise instantons |p5[|). Thus the real question to ask is not why we write the instanton in 
this form, but rather - what is the physical meaning of these parameters, and what is their relation to the original 
parameters characterizing Ho and Hint^. 

We start by describing the adiabatic contribution to the action due to 5k- This comes from the environmental spins 
strongly coupled to the central spin, i.e., with couplings ojk^ ^o- The slow (on a time scale defined by ri~^) rotation 
of the central spin can not cause any transitions from the ground state for these spins (in quantum tunneling problems 
we deal with temperatures T <C f^o - otherwise the " classical" over- barrier transitions will prevail) , and they have to 
follow the direction of the local field acting on them. It is appropriate to include these spins in the " combined central 
spin" . In general, we increase the moment of inertia which will result in a correction to the instanton action. Let us 
calculate this correction using the particular form of Ho{S) given by ( ^.l|) . The corresponding Lagrangian written in 
terms of angles (0, Lp) has the form 

c{e,ip)^is'ipsined + u{e,ip) (3.15) 



?7(6I,(^) = i^ll[sin2 6'(l-f A^sinV)-2/isin6lcos^] , (3.16) 

where = K^/K\\ is the anisotropy parameter, and S" — 5 = 1/2 • Nad gives the number of environmental spins 
locked in a joint rotation with the central spin. For future discussion we also include a weak magnetic field acting in 
the x-direction h = ^^HxS / {2K\\^) ^ 1. As usual P,pOt, we perform integration over the small deviations of the angle 
Lp around zero, and obtain the effective Lagrangian for 9 as 
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The classical equation of motion starting at 9 = 7r/2 at t = 0, and approaching the angle arcsin(/i) as i — > oo has 
the form 

.^^^4(.-sin.)(l-^)^/^ (3.18) 

where fio is the bare bounce frequency (^^). For small h and a large value of the anisotropy parameter (the action 
is proportional to 28/ X]/^, so that very large values of are necessary to make any discussion of mesoscopic spin 
tunneling physically meaningful) we may drop the last factor in the r.h.s. of equation of motion. The solution is now 

„ , , 1 + h cosh riT 

cosh^^. + /. ^'-''^ 

n^^o^Vl- ■ (3.20) 
Substituting this expression back into the effective action we have 



A = 2^JKi\/KjlS'{\/i - /i2 - /iarccos(/i)) . (3.21) 
From this we derive the correction to the instanton action due to the increase of the moment of inertia 



SA,n^^K\\/K^Nad, (3.22) 

which is proportional to the number of spin-1/2 "satellites" added to the "bare" central spin. The result is quite 
general, although the coefficient in front of Nad will depend of course on the particular form of Ho{S). 

For the hyperfine interaction we do not expect any variation in the potential energy. In a more general case, i.e., 
for the dipolar coupling between the S and {(fk}, or when there are interactions acting on {ak} other then their 
interaction with S, the energy U{6, (p) will be modified. This problem is hard to deal with exactly because now the 
potential energy will depend on the orientations ^y/ry between the environmental spins and those forming the central 
spin. However to estimate the effect we do not need the exact solution. It is sufficient to look at the correction to 
the action due to weak applied field, which also changes the potential energy, by amount SU ~ jeHS. It follows from 



(3.21) that for small h the correction is 



SU 

SApot TT-^ , (3.23) 

and may be quite strong if 6U is large. Still we assume that only a small number of " loose spins" on the surface of 
the grain Nad ^ S have dipolar couplings strong enough to move adiabatically with the grain magnetisation; thus 
the total change SApot is much less then Aq. 

To summarize, the adiabatic parameter Sk can vary from a very small value ~ ^^^^^ to ttSU/Qo depending on 
the model ( for nuclear spins inside the grain the former case is the most relevant one since the dipolar coupling is 
much less than the contact hyperfine interaction). 



Now we proceed to the discussion of the terms given by Pk (3^). Suppose we have some environmental 

spins which are coupled to the central spin weakly, Uo/^o <C 1. These spins can not follow the central spin adiabatically 
- they rather see the rotating S' as a sudden perturbation, so that most of them simply conserve their spin orientations 
during the time On the other hand their coupling to the central spin slightly changes the actual instanton path. 

We may start by considering this influence as a weak static magnetic field due to the {ak} acting on 5^, i.e., we write 

the interaction Hamiltonian as Hint = S ■ 6H, where SH is an operator in the environmental spin subspace, and can 



be easily found for the hyperfine interaction (2.4) and dipolar interaction (2.5). Using the fact that S is rotating so 
fast that SH hardly changes during the transition time, we may neglect its time dependence (which will be accounted 
for later in our discussion of the parameters ak and 4>k)- The rest is simple now, for we have to find a linear in SH 
correction to the action, i.e., the variation of A due to the extra potential energy Sa ■ SHa- In fact this was already 
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done for the Ho{S) described by ( |2.lD , and SH afong the af-dhectfon. The linear in h correction is easy to obtain from 



where is the contact hyperfine interaction. In this particular case we find 

6 = --^; Vk^x. (3.25) 

One might try to argue that the calculation of the first order correction to A can be done for arbitrary Ho{S) and 
Hint{{dk\,S) as follows. In the expression for the instanton action 

A = - J driCoir) + H.ntir)) , (3.26) 

the first order correction is given simply by 

5A = - [ dTH^uiSoir)) , (3.27) 



where So{t) is evolving according to the unperturbed semi-classical solution. For the case of hyperfine interactions 
considered above this expression takes the form 

SA = -^^k [ drS^ir) , (3.28) 



25* 

Substituting here ( 3.19| ) (for ft, = 0) we reproduce the result ( [3.25 ). The symmetry of the problem in zero external 



field implies that that the corrections due to SH^.y are zero, which means that here /J^. w 0. 

This would be a correct way of calculating things if we were dealing with ordinary particles but it is not correct for 



the angular momenta. The formula ( 3.27 ) takes into account the change in the potential energy only, while for the 
constrained rotation of S any time-reversal symmetry breaking field can also influence the topological phase. To see 
how it works we apply the magnetic field along y-direction, which results in a term 

- 2is:||/isin6'sin(p w -2/'S:||ftsin6' . (3.29) 

Integrating over small ip as before we obtain the effective Lagrangian for the angle as (up to a constant term) 

c{9) ^ Coie) ^ t^e , (3.30) 

which results in a topological phase correction 

iirS — > iirS + in ^'^ ^ — = inS + iTTT^^^y ■ (3.31) 



Thus we find one of the contributions to the coupling ak to be 

(top) _ 7I"Wo , ^{top) 



y, (3.32) 



Now we turn to the last contribution to our instanton operator (p.6|), which is due to the small, but unavoidable 
dynamics of {ct/c} during the transition. To define this we note that from the environmental spin point of view the 
rotating central spin is seen as a time dependent external magnetic field. Thus we calculate the evolution of the 
microscopic spin wave-function according to 

fk\Xenv)^e^'^"-'^^Uxenv) , (3.33) 

For the weak coupling case ujk/^o ^ 1 we are considering here (sudden perturbation), this can be calculated as 
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fk~l + i5H j drSoir) « exp{i(5i? j drSoir)} , (3.34) 

Thus we find exactly the same expression as in ( 3.27| ) which is hardly surprising since we are dealing with the same 
Schrodinger equation but now in the real time domain. From this we derive the second contribution to the parameter 
ttfc. Again, for contact interactions between the spins wc have 

Sfc = ^ ; nk=x , (3.35) 

These results for ak, £,k, and Sk are all for the realistic case where cok/^o ^ 1- However in the rare cases where 
Wfe/rio > 1, the generalization of the above method is easy. For example, in Fig. 4 we show the variation of ak, as a 
function of cok/^o, for the Ho{S) of ( ^.l[) . The actual form of a{uJo/i^o) for a contact hyperfine interaction is quite 
lengthy, but almost indistinguishable from the simple form 

a{x) - (7r/2) tanhx . (3.36) 

Remarkably, this simple form is almost universal - it hardly depends on the detailed form of Ho{S). This point is 
dramatically illustrated if we take the really pathological instanton 

0{r) = ^ivir + 7T/no) ~ v{r " ^f^o)] + |[r?(r + tt/I],,) + tj{t - ir/no)] . (3.37) 



where tj^t) is the unit step function. Equation (3.37) describes a "jerk start" and "jerk stop" instanton, with the 
polar angle 6{t) varying linearly in between; thus 

'^os" = , 2i^o2U/2 sink/2(l + ^o/f^o)'/'] , (3.38) 

which is almost indistinguishable from Fig. 4 except for some small "wiggles" around a w 7r/2, when lUq > (these 
arise from the infinite acceleration and deceleration at r = ±T:/flo respectively). 

The saturation of a{x) at 7r/2 when x > 1 is easy to understand - it corresponds to the adiabatic (Berry phase) 
limit, when ak rotates adiabatically with S. Similar calculations may be made of ^fe, (j)k, Sk, etc., for the full range 
< a; < oo of couplings. 

This concludes our discussion of the instanton operator. Each term in the exponent (3.6) has its own physical 
meaning and may be calculated from first principles once the form of Ho{S) and the most relevant interactions are 
identified. Here we analyzed the most important case of hyperfine interactions between the central spin and nuclear 
spins when Ho{S) is given by (|]^). We see that the expression ( |3.6| ) simplifies a lot for nonadiabatic coupling 
(iUk/no < 1): 

4 = /3fc = 0fc = ; akfik = ^^(^' y) ; ^kVk ^-^^ ■ (3-39) 

In the adiabatic limit {ujk/^o 3> 1) we have 

;9fe - 0fc = 6 = ; Sk = \~/'^ ; ak ^ 7r/2 . (3.40) 



IV. THREE IMPORTANT LIMITING CASES 



The effective Hamiltonian i?e// in ( 3.10 ) has a very rich behaviour in general. To appreciate this behaviour it is 
best to lead up to it by first studying some limiting cases, each one of which brings out an important aspect of the 
physics. This is done as follows. First we assume that wj! and are zero - although this is a rather artificial way of 
doing things, it enables us to see "pure topological decoherence" , without interference from anything else Next 
we suppress topological decoherence, by making ak very small, and we also make all the ujk equal; however the ratio 
iO^/u)\ now becomes the crucial parameter, governing the behaviour of "pure orthogonality blocking". Finally we 
study " pure degeneracy blocking" ; this is done by making — s- (suppressing topological decoherence) , and letting 
Ljj: = (which suppresses orthogonality blocking) but now allowing a distribution of different values of Uk ~ wj! . We 
shall see the precise meaning of these appellations as we study each limiting case. 
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Naturally, the generic behaviour of Hej / , with completely arbitrary values of ak , , i^l and involves all three 
mechanisms; it can also involve nuclear spin diffusion. We shall see how this works in section 

It is crucial in what follows ro remember that a necessary (but not sufficient) condition for coherence is that the 
total effective bias acting on S be very small - the initial and final states of S must be within ~ Aq of each other, 
otherwise tunneling cannot occur at all. Even if it can, it may be incoherent. 

A. Topological decoherence 

Formally the case of pure topological decoherence applies to the following special case of He/ / : 

N 

HlYf = 2AoTa= COS [$ + ^ akUk ■ <Jk] , (4.1) 



k=l 



Notice here that (a) since tj]! and uj-^ are zero, all environmental states are degenerate, and so (b) the initial and final 
states of the environmental spins, and S, are degenerate - there is no exchange of any energy between S and the {^k}- 
The only thing that is exchanged is phase, i.e., the phase of S becomes entangled with that of the {ak}, during the 
tunneling of S. Thus the initial and final states of the spin environment are not the same. The term proportional 
to ^k which changes the real part of the instanton action and tends to choose either clockwise or counter-clockwise 
trajectories will be included in section 

The properties of -ff*// were entirely elucidated in refs. ||l|,|| (see section 2.2 (ii) of ref. here we explain the 

derivation and results in more detail. What we are interested in are the correlation functions of S{t), most particularly 
in Pf^^(t) and its Fourier transform in the frequency domain. This function simply tells us the probability to find the 
system in the same state at time t as it was at time t ~ 0; its Fourier transform tells us about the response of the 
system to a continuous perturbation at some frequency lo. In the case where all interactions are zero we simply revert 
to the free central spin, for which 

Pu (t) = ^ [1 + cos(4Aoi cos $)] (4.2) 



X'^iij) - nSiiv- I 4AoC0S$ |) ; {for lu > 0) 



(4.3) 



Here x" is the imaginary part of xi^), the Fourier transform of P^^{t). In this non-interacting case $ = nS, and 
Xo is just a sharp line at the tunneling splitting energy | 4AoCOS7rS' |. The spin oscillates coherently between | S^i) 
and I S'2). The series (4.2) is generated from ( [l.l| ) (or alternatively, from the bare transition amplitude (|3.5D), in the 
absence of the {ct/c}, by the usual summation of all possible instanton flips in a time t (see refs. [|l|,|2[33 34[43t ) . 
Returning now to H^^ffi that we can formally write the solution as 



Pni^) = 2 r ""^ [4Aotcos ($ - 



N 

E 

fc=i 



(4.4) 



simply treating the cosine function in as a c-number, because there are no other terms in the Hamiltonian, and 

it obviously commutes with itself. Or, if one prefers an expansion over all instanton trajectories, we can write 



Pi 



n 




\2s 



2s 



{2s) 



^ ' (2s)! ^(2s-?i)!n! 

s=0 ^ ' n=0 ^ ' 



TV 



k=l 



(4.5) 



(...) is a thermal average over environmental states. Here all these states are degenerate, the average is independent 
of T and one finds 



N 



F{n~m) = (Jl 



,iaknk-Sk(2s-2n)\ 



k=l 
N 



cos[2Q!fc(ri — to)] 



(4.6) 



k=l 
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In the non-interacting case, i^(j^) = 1, and we get (4.2). For small ak we may approximate the product in (4.6) as 



N 



(4.7) 



fc=i 



Notice that A is just the mean number of environmental spins that are flipped. We now employ a very useful 
identity 



{2sy. 



to rewrite the expansion (|4.5|) as a weighted integration over topological phase: 



Pn= E ^A(m) rge^^™(*--)ji + icos(2A<,(^)t)| 

m— — oo Ok J 

y m— — OO ) 



(4.8) 

(4.9) 
(4.10) 



where Ao{(p) — 2AoCOS(/3, and J2m{z) is the Bessel function of m-th order. For F{m) — 1 this is just another 
representation of the coherent dynamics ( |4.2| ). On the other hand in the adiabatic strong-coupling limit, where all of 
the 7r/2, we simply have F{m) — (—1)™, so that 



Pn (0 ^ [1 + cos(4Aoi cos $)] ; {at -> tt/2) 



where the renormalized phase is just <& -I- Nn/2, that is 



N 



k=l 



where (j)^ is the Berry phase, defined as before via 



hk - T^/2 . 



(4.11) 



(4.12) 



(4.13) 



Equation ( 4.11 ) simply tells us that in the strong coupling limit the total topological phase now comes from the 
bare phase ttS*, plus the adiabatic Berry phase accumulated by all the environmental spins rotating adiabatically with 
it. If the coupling Hamiltonian is such that, e.g., all the {ct/c} are forced to be parallel to S in this limit, then this 
extra phase will just be N-k/2 (i.e., we have effectively increased our spin quantum number from S" to S" -f N/2). 
However this is by no means necessary - if there is some other field acting on ak in addition to that due to S then 
the Berry phase will deviate from tt/2. 

By far the most important regime is the "intermediate coupling" regime, in which ak > N^^^^ (note that if N is 
large, then ak may be very small here). This is typically going to be the case, since usually ak will be independent 
of N. In this case we have a large parameter A in (4.7) with 



Fxiv) — S^fi + small corrections (intermediate) 



(4.14) 



so that, very surprisingly, we get a universal form in the intermediate coupling limit for P^^(t) (again, ri(x) is the 
step function) 



Pff{t) 



l + JoiAAot) 



(compare the angular average of the coherent series in ( 4.1C| )), and x"i^) 

2 



X"(c^) 



(16A2-w2)l/2 



r](AAo — uj) (intermediate) 



(4.15) 



(4.16) 
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We plot these universal forms in Fig. 5. The physics of this universal form is simply one of phase cancellation in the 
expression (4.6) for F{m). As explained in ref. [Q, this phase cancellation arises because successive flips of S cause, in 
general, a different topological phase to be accumulated by the spin environment, so that when we sum over successive 
instantons for 5*, we get phase randomisation and hence loss of coherence. This forces 2s — n = n in the sum in ( [f.5| ), 
i.e., the only paths that can contribute to P-\^-\^{t) are those for which the number of clockwise and anticlockwise flips of 
S are equal - in this case the topological phase "eaten up" by the environmental spins is zero. By looking at expression 
( 4.15 ) one sees another explanation - the universal behaviour comes from complete phase phase randomisation [Q, 
so that all possible phases contribute equally to the answer! The final form shows decaying oscillations , with an 
envelope ~ at long times. This decay can also be understood M by noting that the "zero phase" trajectories 

that contribute to now constitute a fraction (2s)!/(2'*s!)^ ~ s~^'^ of the total number of possible trajectories, 
where s ~ Kot. Because of this decay, the peak in the spectral function &i oj ~ 4Ao cos $ is now transformed to the 
spectral function of a 1-dimensional tight-binding model, as shown in Fig. 5b. 

From the point of view of fundamental principles this result is really rather interesting, for here we have a case of 
" decoherence without dissipation" ; the coherent oscillations of S{t) are destroyed without any exchange of energy with 
the environment. This mechanism would persist even if the bath were at zero temperature. Whilst there is no reason 
why such "pure decoherence" cannot exist elsewhere in nature, we are not aware of other examples |44 . 

It is actually very interesting to examine the intermediate coupling regime a little bit more closely, especially for the 
1/2-integer spin, when there is no tunneling for a/j = 0. In Fig. 6 we see that weak coupling (uk taking random values 
between and 0.05) to 40 spins completely changes this behaviour - we now have an effective splitting Ae// ~ Ao/3, 
but with strong damping . Fig. 7 shows rather stronger coupling (a^ having random values between and 1). The 
result shown is for S =integer, but in fact the result for 5* = n + 1/2 is almost indistinguishable from this. We see 
that the result is, within numerical accuracy, completely indistinguishable from the universal behaviour in (4.15), for 
only 20 environmental spins. 



B. Orthogonality blocking 

Orthogonality blocking arises when there is a mismatch between the initial and final state wave-functions of the 
environment iPi. Expressed in this way it is an old idea, as old as the Debye- Waller or Franck-Condon, or polaronic 
factors for superohmic environments, or the "orthogonality catastrophe" of the X-ray edge and Kondo problems [^ , 
which reappears in the discussion of Ohmic environments in the oscillator bath model of the environment [p^ . 

In the present case of the spin bath, the mismatch between | a™) and | u^™) arises because the initial and final 
fields, 7^^'' and acting on au, are not exactly parallel or antiparallel. Thus if 3k is initially aligned in 7^^\ it 
suddenly finds itself, after S flips, to be misaligned with 7^ '. There is then an overlap with a flipped state in the 
new basis, equivalent to an amplitude (ik for to have actually flipped (semiclassically, | cr^*") starts precessing in 
7^ ). Depending on the readers' taste, one may think of this as a "quasi-flip" generated by the basis change, or a 
real flip; however quantum-mechanically, for all purposes the spin has actually flipped. 

Thus superficially we are dealing with a standard orthogonality problem. However as we now see the details are 
rather different. 

We start by applying the following restrictions to -ffe// in ( 3.10| ): 



(i) We assume a\N <C 1, i.e., extremely small ak', this removes topological decoherence. 

(ii) We assume that the effective fields 'y^^ and acting on ak at the beginning and end of the instanton are not 
either parallel or antiparallel; usually we will assume however that they are nearly parallel or antiparallel, i.e., that 
either 7^^^ « 7^^^ (i.e. that uj-^ ^ uj\), or else ^^^^ ~ —^^^ (so that cj^ <C wjl). Both of these cases are very common, 
indeed almost the rule - quite apart from any stray fields that might disrupt perfect parallelism, there is also the 
contribution coming from barrier fluctuations (cf. Appendix A) which will contribute to orthogonality blocking even 
if parallelism is otherwise is perfect. Stray fields will come from a variety of sources (weak external fields, or strain 
fields in the sample, or demagnetisation fields from the dipole forces, etc., etc.). 

(iii) We assume that the couplings Uk satisfy the inequalities Wk S> Ae//, where Ag// will be defined below. 
The parameter Ag// < Aq in general; often Ae// <C Aq. This condition is virtually always obeyed for hyperfine 
interactions, since we expect Ag// < IMHz. 

(iv) We assume that all the ujk are equal. This is of course a very unrealistic assumption (even the very small spread 
in hyperfine couplings will turn out to be important), but it is made to avoid the introduction of degeneracy blocking, 
to which we will come presently. 
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To be specific let us consider the case where 7^^^ ~ ^l^k^ which will arise in many situations; the analysis for the 
ther most common problem, where 7^^^ « 7] 
blocking" applies to the effective Hamiltonian 



other most common problem, where 7^^^ « lk'^^ goes through in exactly the same way. Thus "pure orthogonality 



Heff = 2K^fx + Tz-^ ^ ^fc ■ CTfc + ^ TOfc • CTfc . (4-17) 

k=l k=l 

A<i,=AoCos$, (4.18) 



where lOo ^ wj^; in (4.17) we use a "Kramers renormalised" splitting A$. 

We now make a number of observations, which are crucial to understanding of orthogonality blocking: 

(i) The spin bath spectrum is split by wjj into "polarization groups" of degenerate lines, having polarisation AA^ = 
TVT _ 7v^, and energy Et-{AN) = r2w||AA^/2; that is, we have + 1 polarisation groups each separated by energy 
Jl, whose energy changes (by J), AN) if S flips: E^{AN) = -Ei^{AN) = cjJjzW/2. There are C'(^+^)/2 degenerate 
states in polarisation group AA^. 

(ii) For S to flip at all, we require near resonance between initial and final states, i.e., they must be within ~ A$ in 
energy. Since Uo ^ A$, this means that if no nuclear spins are flipped during the transition of 5, only states having 
ZW = (the "zero polarisation" states") can flip at aU (since E-^^iAN = 0) = ^^.(ZW = 0) = 0). At high T, only a 
fraction / — ^J^jixN of grains in an ensemble will have AN — 0. 

(iii) However the presence of uj-^ means that some spins are flipped (in general a different number of them in each 
transition of S\ Now consider a grain in polarisation state ZW — M. If when S tunnels, at the same time nuclear 
spins also flip so that AA^ — ZW, then since £"^(1/) = E\y{—M\ resonance is still preserved and the transition is 
possible. Notice however that (a) for this change in polarisation of 2Af, at least M spins must flip, and (b) if we wish 
to preserve resonance (i.e., S makes back-and- forth transitions) then the polarisation state must change by ±2M each 
time. In fact S cannot flip at all unless there is a change in polarisation state of magnitude 2M. 

These 3 observations will lead to a fourth, which is most easily given once we have defined a function PM(,t), which 
is simply the correlation function P^-f^{t), now restricted only to those systems in an ensemble for which ZW = M. 
Since, as remarked above, we can only have transitions between ±ZW states, we can write the complete correlation 



function for an ensemble of grains described by (4.17) as 



N 



Puit;T)^ J2 w{T,M)PM{t), (4.19) 
w{T, M) = ^-i^f +*^)/2e-A^-i/'=-i^ , (4.20) 

where Z is the partition function. 
We now make the observation 

(iv) If we are interested in the incoherent relaxation of particularly at long times, then obviously all the Phiit) 



in (4.19) must be considered. However for coherence it is obvious that the Po{t) must completely dominate, since it 
involves transitions of S with possibly no flipping of the {f?fc}, which are necessarily coherent. On the other hand 
PM^o{t) not only involves many spin flips each time (necessarily causing loss of coherence), but it also involves only 
processes in which some number > M of spin flips leads to exactly a change 2M in polarisation. From this argument 
we see that PM^oit) will describe slow (i.e., for large M on time scales ^ Aj^) and incoherent motions of 5*. 

These observations are borne out by the calculation of PM(t) in Appendix A. 2. Assuming that luI 3> lu^ as above, 
one easily sees that the flip amplitude f3k is just the small angle defined as shown in Fig. 8, via 

cos2/3fe = -7^i).7f , (4.21) 
where 7^*^ are unit vectors; it is then found that 



Pm(0 ^ dcx e^'^' ^1 + cos[4A<i, JAf(2\/^x)i]^ (4.22) 



00 



2/ ckxe-^'pll\t,AM{x)) , (4.23) 
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where P^^ {t, Am) is just the free function (Eg. ([4. 2D), but now a function of an x-dependent sphtting Am{x) 



tt 

2A$ Jj\/(2y^x); and k is just the orhtogonahty exponent, i.e.. 



N 



e-'' ^llcoaPk , (4.24) 



fc=i 



so that K sa l/^X^fe-^fc- Thus the frequency scale of Puit) is set by Am (a;). The terms Puit) are easily verified to 
be incoherent, and so 

Pni^) ~ fPo{t) + incoherent , (4.25) 



where / = ^j2/nN as before; the exact answer is given still by (4.19) above 



The frequency scale of oscillations in P{){t) is easily found; we write it in terms of a renormalised splitting Agjy: 

Po(i)-l/2[l + cos(4A*e-«t)] 1 ^ .4 261 

Ae// = 2A$e-« J ^ • ^ 

for small k, and 



Po{t) = 1 - 4(A|tV(7r«)i/2) + 0{K%t^) 
Ae// -2A$/(^k)V4 



(4.27) 



for large n. Thus the naive "orthogonality catastrophe polaronic band narrowing" argument to derive Ae// only 
works for k <C 1 (where it is superfluous!). In the interesting regime k ^ 1, we get a quite different answer (with 
a much smaller suppression). This is understood as follows. In the usual oscillator bath models, band narrowing 
comes essentially without any bath transitions (most of the polaron "cloud" is in virtual high-frequency modes) - it 
is adiabatic. Here, however, roughly k spins flip each time S flips (the probability of r flips is n'^e^'^ /r\, which peaks 
at r ~ At), even though we only consider Po(0: i-^-, even though AN — and does not change (the flips compensate - 
there are just as many one way as the other). 

This inevitable flipping of the {ak} means that even Po{t) is incoherent if k ^ 1. This is most easily seen by looking 
at Xo(^) which has extraordinary behaviour. One finds 



XM=o(^)-\/-r7-T^E 



T^N 4A<i.y^^ |Ji(2V^a;j)| 



(4.28) 



The first factor is just the statistical fraction of zero polarization states in the spin environment. It is the complicated 
structure of the set of roots Xj to the equation Jo(2y^Xj) = ±w/4A$ that yields the interesting behaviour. For small 
K only the lowest square root Bessel function singularity appears, coming from xi ~ [(1 — ti'/4A$)/K]^/^. This leads 
to 

Xm=o(^) = cxp I '^ "/^'^ ) r?(4A$ - w) , {k « 1) , (4.29) 

4A$K y 4A$K J 

and almost coherent dynamics for S*, as we see in Fig. 9. However as k increases, more and more square root singularities 



contribute to x" (which also gets pushed to lower w, according to (4.27)), and it develops the bizarre structure seen 
in Figs. 9& 10. This structure has quite counter-intuitive consequences for Po{t)] interference between the various 
peak structures can conspire to give the correlation function for S a multi-periodic behaviour. An example is shown 
in Fig.ll. 

Readers who wish to more fully understand the mechanism of orthogonality blocking are urged to study Appendix 



A. Recall also that for a real system, the values of the Pk will be renormahzed by the /?^, (Eq.(3.10) et seq.), which 
describe the effect of barrier fluctuations. Thus the real value of k, and the real values of uoj: and a;|, will already 
include these barrier fluctuation effects. 
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C. Degeneracy blocking 



The previous 2 limiting cases studied have imposed some rather artificial restrictions in order to bring out 2 physical 
mechanisms involved in decoherence from spin environments. We now study a third limiting case which introduces a 
note of realism. To do this we suppress topological decoherence (by letting <C 1) and we suppress orthogonality 
blocking (by letting = 0) . On the other hand we now allow a spread of coupling energies wji around a mean value 
Uo- Thus we consider the effective Hamiltonian 

N 

Hef f = 2A*f, + f, ^ a| ; (4.30) 

with a spread of values of c^l of 



y/^(4^ - ^o)' = N'^^SiJo . (4.31) 

i.e., a distribution of width ScUq- Clearly, this Hamiltonian is identical to the standard biased two-level system with 
the bias energy e depending on the particular environmental state; thus e = J2k=i^k'^k- introduction of this 
spread is to destroy the exact degeneracy between states in the same polarisation group. For coherence to take place, 
we require the initial and final states to be within roughly Ag^/ of each other, otherwise Landau-Zener suppression 
of the transition takes place (cf. Appendix A.l). Since we assume a random distribution of couplings around lOq, it 
is clear that there will be a crossover between unblocked behaviour for Swo <C Ao/N^^'^, and Landau-Zener blocked 
behaviour for SlOq ^ Ao/N-^^^. In practice once SlOq ~ Ao/iV^/^, the crossover is essentially complete, and tunneling 
is blocked for almost all grains in an ensemble. Thus for large N, degeneracy blocking is a very powerful mechanism 
for suppressing coherent oscillations of S, since only a tiny spread SlOq will be required to do the job. Now we look at 
this quantitatively. 

The trivial case is of course where SlUq = 0, and the spectrum of {tTfc} is organized in highly degenerate lines 
corresponding to different polarizations e = lOqAN /2. Then, as we discussed in the last subsection any zero polarization 
state of {fffc} gives an exact resonance for S*, and the corresponding correlation function is an ideal one , i.e., we get 



P'f^^it) as in (4.2); however only a small fraction of possible environmental states have AA^ = 0, and for large N this 



fraction is defined by the central limit theorem as / = ^J2/ttN . Thus the corresponding spectral function is like ( |4.3| ) 
but with a total weight /, i.e.. 



^^{uj) ^ ^2^/N5{uj - 4A$) ; (5w„ - 0) . (4.32) 
We may simply disregard the contribution to x" coming from states with nonzero AN , because for ujo ^ ^e// these 



states are too far off the resonance. Thus, as we see from (4.32), in this special case the only effect of the coupling 
uJo to the N spins is a reduction of the spectral weight - otherwise the central spin behaviour is completely coherent. 
This is hardly surprising, since we have removed all decoherence mechanisms - we have no topological decoherence, 
no degeneracy blocking, and no orthogonality blocking! This case was also studied very recently by Garg Q in an 
attempt to discuss experiments in magnetic grains. 

Going now to the case where we do have some degeneracy blocking, it is useful to first study the regime where 
/i = 6lOoN^^'^ /u>o > 1 (in real situations, this regime will almost always prevail in any case). Under this regime the 
different groups of environmental states will be spread out and overlap, and we will then get a continuous Gaussian 
distribution for the bias energy e: 

= (0 \ni/2 exp{-eV(2c^o^)} ; (4-33) 
tt'o(27rA')^/^ 

N 

e^Y.4'^k- (4.34) 

fc=i 

This situation is depicted schematically in Figs. 12a and 12b; Fig. 12b also depicts a case where /i < 1, i.e., with 
very small degeneracy blocking. The correlation function P^^(t) for an ensemble of central spins is now obtained by 
averaging the known analytic behaviour for the simple biased 2-level system. At high temperatures 
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Puit) = / deWie 



2A| 



1 - cos(2<y^^74A|)) 



(4.35) 



oo 

1-2A^ J2fc+i(4A$t) . (4.36) 



A = 27rA<i,/(2woV2^) = 27r/ • A* /AiVo ■ (4.37) 



and a spectral function 



8 A 

X"i^) = ^ ; ^(^ - 4A<,) , (4.38) 

;2 - 16A2 

which is shown in Fig. 13. (Integrating over e we used the inequality W^(0)A$ <C 1). We see that not only is there 
line-broadening, but also the amplitude of the oscillations, for the ensemble of central spins, is further reduced by a 



factor A^/2lOo, as compared with the trivial case in (4.32), and by a total factor A as compared with the ideal free 



central spin in ( |4.2| ). We should note at this point that in almost any realistic situation, a random distribution of 
values for w| for a single central spin will lead to results identical to the above, i.e., the ensemble average here is 
merely a device to define the distribution W{e). The only difference is that for a single central spin one has to consider 



(4.36) as a statistical average over many experimental runs during a time which is longer then the NMR relaxation 
times in the nuclear subsystem. 

Turning now to the general case, where the parameter n = SuJoN^^^ /ujo can take arbitrary values, we notice that 
the integral over e in ( 4.36| ) is convergent on a scale defined by A$, and this result will not change at all for any 



distribution which is flat around zero. Thus for lUo 3> ScJqN^^^ > A$, i.e., <C 1 we have the same answer as before, 
but the amplitude A is now replaced by ^ = / • A^/{2\/2ttN6u!o) — A,s,/ {2ttN6uJo)- Thus we find that even a very 
small spread in the coupling constants will be sufficient to broaden the spectral line. Note that neither the energy 
dissipation nor the phase randomisation are responsible for this behaviour, which is totally due to the statistical 
average over the initial states of {^k}- 

V. THE GENERIC CASE 

As we shall see, if we combine all 3 mechanisms just studied, and calculate PM{t) for the full H^ff, we get a very 
complex (but still a closed form) answer. Thus in order to make the derivations more transparent, we first combine 
the mechanisms in different combinations. This tactic is not just pedagogically useful - it also gives 3 more limiting 
cases which correspond to experimentally realistic situations. All this is done in Se ctions 0.A - 0D. Finally, in^E, 



we deal with the residual inter- nuclear couplings mentioned in Section III, Eq.(3.12). Within a particular polarisation 
group AA'^, the inter- nuclear dipolar couplings cause "diffusion in bias", simply by flipping random pairs of nuclei. We 
study the role this has on the general solution for the coherent part of Pff{t) given in 

A. Projected Topological decoherence 



We start from the case of "projected topological decoherence" , defined by the following special case of (3.1C) 

JV N 



Heff = 2Aof^ cos [4' + ^ akfik ■ ^k] + ujk^k > (5-1) 
fc=i fc=i 

with the restrictions: 

(i) The effective fields 7^^^ = -7^ ^ , so cj^ = 0. 

(ii) The couplings Wk 2> Aq , and all couplings are equal (thereby eliminating degeneracy blocking). 

The second condition forces the polarization of the environment with respect to the central spin orientation to 



remain constant. Now at first glance we have no orthogonality blocking at all in (5.1), since = (compare section 



IV.B), and so it looks as though only zero polarisation group states will have any dynamics. However this is incorrect; 
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because now the dynamic couplings ak can flip spins when S flips, with amplitude ak] they play exactly the same 
role as the f3k in pure orthogonality blocking (Section IV. B), and so we can say that in ( ^.l[ ), orthogonality blocking 
is generated by the ak- Thus, in the same way as in Eqs.( 1.19 ) and ( 4.20 ) we write Pffit) as a weighted sum over 
PM(t), where now (cf. Appendix B.l): 



PM{t)^ I dxxe-''' Pyi^) f — e*^"^*"^^ {l + cos[4AoUM(2a;%/A^)cos(p]} 

= fdxxe-^'ll+ ^A'(™)e*""*J2m[4Aoi 

m— — oc 



,tJM{2xVX^)] 



where F\i{m) is defined analogously to F\{m) in (4/7), viz. 



(5.2) 
(5.3) 



Fy{m) = e 



-4A'r 



(5.4) 



1 ^ 



(5.5) 



fe=i 



and we notice that A > A'. Eq.(5.2) can be interpreted in 2 ways, both of which are important. First, we can interpret 
it as an orthogonality-blocked expression, now with frequency scale 



Am{(P,x) = 2AoCos(i^)Jm(2x\/A - A') , 



(5.6) 



which is then averaged over ip, to give phase randomisation (compare Eq.(4.9)). Alternatively we can regard it 



is an integration J dx over an already topologically decohered function p||°^^(f, Aj\/(a;)), with A]\j(x) now given by 



Am{x) = AoJM(2a;\/A — A'); compare again with Eq.( 4.23 ). 

This is the result for Pm {t) ; but now we can go through exactly the same arguments as for pure degeneracy blocking, 
to show that only Pq (t) may behave coherently, with a weight ~ y^z/TrN. Again the exact expression is given by the 
weighted sum ( 4.19| ). 

There are various interesting cases of (5.3) for M — 0. If A = A' (i.e., is parallel to z), then we go back to 
pure topological decoherence - the projection operator then commutes with the cosine operator. On the other hand 
if A' = 0, we have pure orthogonality blocking as stated earlier, and in fact when A' = 0, the parameter A plays the 
role of K in ( 4.22| ). Notice that whereas the case A = A' can only occur accidentally. A' = is quite common - indeed 
it pertains to the model in Eqs. ( |6.l[ ) and ( |6.2| ) which we will examine below in connection with real experiments. 

We really begin to see the analogy with degeneracy blocking when A, A' 1; just as with pure topological 
decoherence, Fy (m) collapses to a Kronecker delta, and we get the universal projected topological decoherence form: 



dxxe 



1 + Jo[4Aoao(2xVA- A')] 



(5.7) 



dxxe 



[16A2 J2(2a;x/A^) - l^2]i/2 



r;(4Ao I Jo(2a;VA^) | -iv) 



(5.8) 



which generalizes the result of (4. If) for pure topological decoherence. Eq.(5.7) should be compared to ( 4.22| ). 

We show in Figs. 14 and 15 some results for x"{^) for selected values of A ~ A'. The results are startling; even a 
very small value of (A — A') significantly washes out pure topological decoherence; but for any large value of A', we 
never get back the pure orthogonality blocking spectrum. 



B. Orthogonality blocking plus degeneracy blocking 



We now look at what might be called "biased orthogonality blocking"; as explained in section ^.C, if we introduce 
a spread StOk about some luq, then in effect we introduce a bias e which destroys the degeneracy of the environmental 
states. Thus we now study the Haniiltonian 
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H. 



: 2A 



N N 



(5.9) 



k=l 



k=l 



with the spread in frequencies ojk = [(^H)^ + {^'t)'^V^'^ defined by Eq.(4.3f ). The effective bias e — X^fcLi ^l^k again 



and now we would hke to take account of the effect of this bias on the orthogonahty blocking. Again, one starts by 
calculating a function Pm referring to polarisation subgroup M, but this time in a bias e. One finds (Appendix B.2): 



PM{t,e) = 2 / dxxe-'' \ 1 - 



Ai,,(x) 



1 - cos [2tE{x) 



te--Vjy(t,AM(a;),e) . 



(5.10) 



E^x) = + Aliix) . (5.11) 

This expression could have been guessed without any calculation - it is just the usual orthogonality blocking average 
J dx, but now for a biased system (instead of, as in pure orthogonality blocking, an unbiased expression - cf. ( 4.22| ). 

To get P^-^it), we must now both integrate over bias and sum over M. Thus the generalisation of our previous 
weighted average ( |4.19D is just 

P(t;T) ^ deWie)—- ^ Puit, e ^ Mcj,) . (5.12) 

M=~N 

However, yet again, we use the usual argument that only Pq may produce coherent response. Assuming again that 
the spread in couplings satisfies fi — N^^^SuJo/i^o > 1 (a trivial modification is necessary for /icJo/Ao ^ 1 but /i < 1, 
as discussed in section ^.C), we then get for M — 

/CO 
dxxe-""' I Jo(2V^a;) | ^ Jafc+i [4A<i> | Jo(2V^x) | t] (5.13) 

/oo 
dxxe-=''2A{x)J2J2k+i[2Ao{x)t] , (5.14) 
k=0 

where the a;-dependent spectral weight is 

A{x) = AJo(2y/^x) . (5.15) 

Yet again, by comparing with the previous results for pure orthogonality blocking and pure degeneracy blocking, we 
could have guessed this (compare Eq.(4.36); we have just done an "orthogonality average" over a "biased averaged" 
expression for the free system with a;-dependent tunneling frequency Ao(a;). 

The expression for Xa/=o('^) follows immediately; we have 

X^'(^) = ^ I dxxe-^\A{x)j-^-^^^v{^ ~ 2 I Ao(x) |) , (5.16) 

Figs. 16 & 17 show some representative plots for this coherent part of x"{i^)- it is in fact almost completely incoherent, 
with total spectral weight 



(5.17) 



{daj/2Tr)x"{uj) ^2A drae"^ | Jo{2y/l^x) \ 
_ 2r(3/4) A 

- ^3/2 ^-1/4 ; 



a result which is very accurate even for k ~ 0.02. Of course, once we include all other polarization sectors M 0, 
the total power increases - in fact it now becomes ~ An^^^ for large k, since ~ k^/^ different polarization sectors 
contribute. However none of this is coherent either. 
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C. Pure topological decoherence plus degeneracy blocking 



Here we discuss the case when the central spin dynamics is governed by pure topological decoherence is a biased 
system. With what we have calculated already the answer is readily obtained as follows. We note that the correlation 



function in Eq.( 4.10 ) is clearly given as a coherent correlation function Pi^{t) (see ( [A3| )) integrated over ip and 



summed over m. Now, if we add the bias, all we need to do is to substitute for the unbiased correlation function with 



the non-interacting biased one from Eq.( AlC ). Thus we have for the biased topological decoherence 

-inc., (5.18) 



Pnit) = Ej>^(-^) I {l - (^1 _ eos [2t^e^ + Ag(^)]~ 

Integrating over the bias with the usual assumption about the distribution function we immediately find 

oo „ J oo 

Pnit) = 1 - E / ^e'2™(*-^)2A(^) E J2fc+i[2Ao(^)t] + inc. , (5.19) 

m=-oo k=Q 

for the correlation function, with Ao((^) = 2A$cos(/7 as before, and A{(p) — A cos ip; and 

;^"(^) ^ M y F^{m) [ %^^rai^-^) -?7(c./4A„- | cos^ |) , (5.20) 

' ' ^ ' } 2t: [(w/4Ao)2-cos2(pli/2 / o i , y , 



m— — OQ 



for the spectral function. 

One might be confused whether this calculation applies to any real system, because on one hand the pure topological 
decoherence implies no (or at least very weak Uk ^ Aq) interaction between the central spin and its environment; 
on the other hand including degeneracy blocking we assume that the coupling is strong enough to move the system 
out of the resonance. How one can reconcile this? Below we demonstrate that this situation is described as a special 
limiting case of our effective Hamiltonian. 

Suppose we have a very weak interaction between the central spin and its environment so that ujk ^ A,,. Since 
Ao is exponentially small as compared to fio we are dealing with the case of extremely small a/c ~ Wfc/r^o <C A.o/flo- 
Imagine now that we have a really macroscopic number of environmental spins contributing to the decoherence so 
that the small value of ak is compensated or even superseded by N so that A ~ Na^ is of order unity or even large. 
One immediately understands that in this case it is impossible for the central spin to tunnel and not to flip some 
environmental spins, but since tok is much smaller than Aq we are still in resonance (if < Ao) and may neglect 

energy dissipation. 

It is quite obvious, however, that for such a large number of microscopic spins both the parameter fi = 
SiUkN^/-^ /ujo S> 1, and the level distribution width N^^^u>o ^ Aq are very large, and we do have to consider tunneling 
in a biased system. The reason why the two mechanisms are combined so easily is that the bias is produced by the 
total effect of all spins, while only an extremely small fraction of them flips when the central spin tunnels. As described 
elsewhere the numbers are such that this special case has direct application to the environmental spin decoherence in 
SQUIDs. For large parameter A Eq.( ^.20 ) simplifies to 



„, , 2 A f dip cos^ G9 , , T , ,s . s 

X" = — / 77 C./4A,- cos f , 5.21 

uj J 2tt [(tj/4Ao)2 - cos2 (^]V2 



which can be expressed in terms of Elliptic functions. 



D. Combining the 3 mechanisms: the generic case 

With the preceding studies we hope that the following basic message has come through, viz., that combining 
combinations of the 3 mechanisms simply boils down to combining 3 weighted averages. These averages are 

(a) A "topological phase average" given by ^ ^a(^) / ^e*2;v(<i.->^) . ^5 22) 



I/— — 00 
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2 

(b) An "orthogonality average" given by 2 / dxxe"'^ ; (5.23) 

Jo 

(c) A "bias average" / deW(e)—- . (5.24) 
J Z[P) 

These averages apply to any polarisation group M ; the final expression for P^^ [t) is then given by averaging the 
quantity 

P{°;(t; AM(^,x);e) = 1 - ^p^sm\EM{v,x)t) , (5.25) 



i.e., the free biased correlator for a splitting A 



M- 



AMi(p,x) = 2AoC0S(^ Jm (2a; VA- A') , (5.26) 

El,{ip,x)^e^ + Al,iip,x) . (5.27) 
Thus, using this prescription, we can immediately write down almost the final generic case; we have, assuming that 

'k 

N N 

Hef f = 2AoT^ COS [$ + ^ a^nfe • ffk] + ^ ^ct^ . (5.28) 



in the effective Hamiltonian both uj-t = and £,k — 



2 

k=l k=l 

the full correlator is given by 



N 

P^^{t;T) = I deW{e)^ ^ PM(t,e - MoJo) ; (5.29) 



Z{(3) 



PM{t-e)^2 r dxxe-^" F'^ f ^e^^. 



1 ^ ^^TT 7 sm GBm(^, a;)i) 

Ejj{ip,x) 



(5.30) 



at any temperature. It should be emphasized that this expression is exact provided the ak are small and ojk ^ A$ 
(which is virtually always the case in nature). 

If we are only interested in possible coherence, then it suffices only to look at M = term in this expression, for 
reasons we have previously discussed. Carrying out the bias average then just gives the obvious generalization of 
Eq.(pl|), i.e., 

Pnit) = 1-2 / dxxe-^"2 V F^(i.) / /e'2-(*-v)A((^,a;) V J2fc+i[2Ao(^,x)i] + inc. ; (5.31) 

•^0 ,^=~oo J k=0 

A{ip,x) = Acosip Joi2xVX - A') ; (5.32) 
and an equally obvious generalization for x"(^) 

X'^iu;) dxxe-^' £ F',i.) j ge^^-^C*".) 

z^— — oo 

A{lp,x)Ao{lp,x) 

[c.2-4A2(^,:e)]i/2 ' ^^-^"^^ 

These results for the effective Hamiltonian in ( |3.1Cl| ) are all we will derive here. We now argue that these are all we 
will ever really need for qualitative understanding. We hope it is clear that if we put the Luj: term back into Heff^ it 
would be equally straightforward to calculate P^^(i), since the role of the Pk in the calculations is analogous to that 
of the afc, but even simpler to deal with. The only term we have still not dealt with is the term ^fcWfc • ct^ which also 
appears in the cosine term in H^ff. In Appendix C.3 we show how this can be included - the final expression is a 
monster, but can easily be put onto a computer for comparison with real systems (as in the next Section). We also 
show plots here of Xoi^) '^^ Fig. 18 and Fig. 19 for a few representative values of A — A'; one should again notice how 
rapidly decoherence sets in, and how large is the power reduction factor. 
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E. Nuclear Spin Diffusion 



Finally we comment on the role of the role of the very weak dipolar/Nakamura-Suhl couplings between the nuclei 
(Eq.(|l|). There is potentially very important decoherence effect that can come from these interactions, because the 
dipolar coupling can cause pairwise flipping (|ti) ^liT)) of the nuclei. Even though this interaction is weak (it is of 
strength T2^, and will be in the range 10^ — 10^ Hz - see Table ||) it is important because it allows both real space 
spin diffusion and diffusion of the system in a bias energy e, within a given polarisation group. For a set of N nuclei 
the pairwise flips occur at a rate ^ -/VTj"^; each flip will change the bias energy by ^ SlOq in energy, and the bias will 
random walk throughout the energy range N^^'^SuJo = of the polarisation group. Now suppose AA^ = M , so that 
transitions of S can only proceed if the time it takes e{t) to diffuse out of the energy window around resonance, of 
width Am, exceeds the "waiting time" Aj^ necessary for the transition to proceed. Since in a time A^^^, roughly 
N/T2AM pairwise flips occur, e{t) will change by roughly 5uJo{N /T2AmY^'^ in this time, and so coherence requires 

Ai,»^(5a;o)2. (5.34) 

J-2 

If the spread 6uJk arises itself from Nakamura-Suhl interactions between the nuclei of the same kind, Tj"^ ^ 5oJk, and 
this criterion becomes Am ^ T^^N^^^- These criteria are very hard to fulfill unless M is small, i.e., for the zero 
polarisation group. One obvious way to avoid this nuclear spin diffusion decoherence mechanism is to choose a system 
like Fe or Ni in which nuclear spin isotopes are rare (cf. Table |). Then T^^ will be very small, and coherence may 
be visible for quite large N (and even larger S, of course). 

In the opposite case, when e{t) fluctuates too rapidly for tunneling, we will simply get incoherent relaxation of 
Pmit), with a correlation time tm = (2A|^//xWo)~^; or, deflning a correlation time tm{x) via 

r-/(a:) = ^^M£), (5.35) 

we would find a completely incoherent form for P-^-^ (t) : 

P^^{t) = 2^w{M;T) I dxxe-''\xp{-t/TM{x)} , (5.36) 

(and obvious generalisation, if necessary, to include integration over ip). 

This concludes our analysis of the central spin model. We should emphasize that whereas in this paper we con- 
centrated on the coherence problem, the general formalism we have constructed, and the results we established, can 
also be used in analysing any of the low-T relaxation properties of S. The generalisation to a finite applied field, for 
example, can be done straightforwardly,starting from our generic expression. 

VI. A REALISTIC EXAMPLE 

Readers (particularly experimentalists) may feel that in spite of these analytic results, there is still a very long path 
to follow from the bare, untruncated Hamiltonian Ho{S), with some set of couplings to nuclear spins, paramagnetic 
spins, etc., to our flnal answers. We therefore, in this section take a specific bare Hamiltonian, viz., 

/f(5,{/fc}) = i7o(5') + ^5]^fc5-4, (6.1) 

k 

Ho{S) = -K\\sl+K^sl , (6.2) 

with a range of couplings around a central hyperfine coupling lOq between S and the nuclear spins 1^. We then 
solve this problem almost completely, by giving the analytic forms for i^^f (t) and x"(<^)- We hope that this will make 
it clear how one may deal with whatever realistic central spin Hamiltonian that Nature may care to throw at us. 
We observe that this model Hamiltonian actually provides a very good description of magnetic grains (compare the 



examples in section [II. B.) provided we ignore: 



(a) The coupling of the grain order parameter to both nuclei and electron spins outside the grain; 
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(b) The coupling to "loose spins" on the grain surface; 



(c) The breaking of the degeneracy between | 5*1) = z \ S) and | S2) — —z \ S) in (|6.2|), caused by strain fields, 
etc., coming from non-uniformity in the shape of the grain, surface defects, etc. - this will lead to a finite "bare" 
orthogonality blocking term, in which /3k now is given by cos 2/3 = —Si ■ 82/8'^. 

The neglect of these terms may, for a real experiment, lead to much greater decoherence than a calculation based 



solely on (|6.l[ ) and (6.2) can give. Thus the results we will give here will describe what one would see in an experiment 
with the minimum possible decoherence. However we shall see that even this minimum is very large - and therefore 
our prognosis for observation of any coherence in magnetic grains is very poor (as already noted in refs. |p]-H|). 



The decoherence that does come into (3.1) and ( |6.2[ ) is as follows: 

(i) There is the inevitable topological decoherence, and also the effect of asymmetric barrier fluctuations. 

(ii) As discussed in section ^A, even if there is no bare orthogonality blocking in (6.2), it will be generated by the 



solution, via projected topological decoherence. 

(iii) finally we note that there will be a spread 5ujk in the couplings {w^} in (6.1). This is an "NMR linewidth" 
effect, and is accessible experimentally. Thus if one was able experimentally to eliminate the effects (a), (b), and (c) 
described above (using, e.g., a perfect grain or magnetic macromolccule, in a matrix of ultrapure *He-ll), then our 



theory would be testable by direct experiment: all parameters in (3.1) and ( |6.2[ ) would be known from the start. 

The spread will in reality arise from various sources. First, there are the Suhl-Nakamura interactions, which can 
give a spread dtok ~ 10 — 1000 kHz\ in most magnetic materials, this spread is at least 100 kHz for the magnetic 
nuclei; some examples are shown in Table ||. Second, there will be a spreading caused by the non-uniform dipolar 
field from the electronic spins which make up the giant spin (the nuclear dipole-dipole interaction will give a spreading 
^ 1 kHz, and can be completely ignored). This spreading will depend on the size and shape of the grain (it is not 
present in a bulk sample , and so NMR linewidth measurements in the bulk will miss this effect). Its size can be 
estimated from simple magnetostatics to be ~ 100 kHz for small grain [S < 10"^) to above 1 MHz for large grains. 
Thus to find Suk in a set of real grains, NMR linewidth measurements should be done on the same grains. 



The solution to the general problem with the Hamiltonian (6J_) and ( |6.2| ) has already been given; the relevant 
effective Hamiltonian is 



N N 



Heff = 2Ao{f+ cos [$ -t- ^{akfik - i£.kVk) ■ ^k] + H.c.} + T^^Jllk ■ Sk ■ (6.3) 

k=l k=l 

in which the parameters are given by 

Ao = Ao 
akrik = -j^y^^v) 

CfcWfe = -^f, (6.4) 

and the unit vector Ik = z. As just discussed, the parameter ojk (and its spread) is accessible experimentally using 
NMR, as is Aq (which is unrenormalised here, for weak interactions) . If we can determine experimentally the values 
of 7^11 and Kj^, we can determine VLo via flpof 

no^'^{K\\K^f'^ . (6.5) 

or, alternatively, fio can be determined directly from microwave absorption experiments (which "rock" the central 
spin back and forth in one of the degenerate wells). 



Equation (6.3) corresponds to the case of projected topological decoherence with no orthogonality blocking, but 
degeneracy blocking coming from the spread 5uJk . This case was completely solved in section ^; these equations are 
valid for all the examples we will discuss, since they assume /icJo/Ao — N^^^SuJk/ Ao > 1. To see what we get, let us 
look at two examples: 

Example A: TbFe^ grains: We take this example as a representative case in which topological decoherence is very 



important; as noted in section [11. B, one has for this system ilo ^ 60 GHz, and the quadrupolar split lines have 
ak 0.08, 0.07, and 0.05 respectively. Now since no experiments to look for coherence have been done on this system, 
we will simply imagine one in which we choose S = 10^, and Aq = 1 MHz (this value of Aq may be optimistic, but 
is not in conflict with equation (|3.2D, given the uncertainty in the ratio of Ki\/K± for TbFe^). We also note that the 
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spread in au caused by the Fe^^ nuclei will give StUk ~ 50 MHz, so that with Nrb = 2.5 x 10^ we have fj, ^ 1. Thus 
we have the following parameters 

a^V2^ = 0.05, 0.07, 0.08 

(wo/Ao)/i ~ 3 X 10^ (TbFes) , (6.6) 

leading to the values 

A - 10 ; (A' = 0) . (6.7) 

The form of x" {^) for this case is shown in Fig. 20 (again only the M = contribution is shown). We see that there is 
complete loss of coherence, and extremely strong suppression of the absorption power {A ^ 10~^ — 10~^). As noted 
previously in Sec. ^.B adding in the M ^ contributions will increase the total spectral weight, but this contribution 
will be even more incoherent and centered around lu = 0. Nuclear spin diffusion effects will of course change this 
somewhat. 

Example B: Ferritin grains: This is a more interesting example, since topological decoherence is unimportant here. 
We shall first analyze the experiments of Awschalom et al. whilst ignoring the interaction between the ferritin 



molecules. As we saw in section [ILB, ferritin is an antiferromagnetic macromolecule with 4500 Fe ions. It also 
has an excess ferromagnetic moment with a somewhat uncertain value (estimates range between 217 and 640/iB). In 
reality, because of this moment, ferritin molecules will interact. There will also be a coupling to nuclear and electronic 
moments outside the ferritin, which we ignore here. 

We are then left with an Fe hyperfine coupling ujf; ~ 64 MHz; there will also be other hypcrfinc couplings to the 
H^ protons in the H2O in ferritin; we also ignore this for the moment (but see our discussion below). We assume a 
Suk ~ 200 kHz arising from a combination of Suhl-Nakamura and non-uniform dipolar interactions - this is probably 
an underestimate, but the quantity has not yet been measured, for ferritin molecules, as far as we know. 

According to Awschalom et al., the anisotropy constants and Kj^ obey {K\\Kj_y^^ ^ 1.72 Tesla, so we shall 
assume a value for ilo of 40 GHz, as a rough estimate Finally, the experiments of Awschalom et al. claim 

that a resonance they see of frequency just under 1 MHz corresponds to coherent tunneling; therefore we assume 



Ao ~ 1 MHz. We thereby arrive at the following numbers, using ( |6.4[ ), (|6.5|), and the definition of fi: 

a = V2^ - 2 X 10~^ 

^ = N^^^ScJk/ujo 2 X IQ-^ 

{ujo/^o)ti - 2 , (6.8) 

where we note that there will be '-^ 100 _Fe^^ nuclei in the molecule, so iV^/^ — 10. Now from (6.8) we find 

A - 4 X 10"^ ; (A' = 0) ; [Ferritin) , (6.9) 

which rules out topological decoherence; since we assume Pk = (per fect sample) this leaves only degeneracy blocking, 
and since {dq/ ^o)lJ' > 1 we can immediately apply the result ( 4.3^ ) in the text, which is depicted in Fig. 13. If one 



compares the curve in Fig. 13 with the data of Awschalom et al. |17|, one is immediately surprised that the lineshape in 
the experiments is not at all that different from the calculation. Now it is important to understand that the reason for 
the sharp theoretical lineshape in Fig. 13 is precisely because topological decoherence and orthogonality blocking have 
been dropped. It will be immediately obvious, on referring to Table I, that the value of the dimensionless topological 
decoherence parameter A is extremely low for ferritin; the value of 10 for ThFe^ is more normal. If we can ignore 
any orthogonality blocking, then it rather seems as though the experimentalists have hit upon an almost ideal system 
for the investigation of MQC in magnets (it could be made even better if one could isotopically purify ||^ the Fe in 
ferritin, to get rid of all Fe^"^ nuclei). Do our calculations then support the thesis of Awschalom et al., to have seen 
MQC in ferritin? At this point we feel that one must be cautious. First, we notice there is a strong reduction, by a 
factor A — Ao/ (2uJoV2ttN) ~ 3.5 x 10"^, in the total spectral weight of this peak, as compared to the situation where 
we ignore the nuclear spins (notice that the question of power absorption in these experiments has already caused 
some controversy ]l8| , p6| ). This absorption is a consequence of the restriction to the M = polarization sector, and 
the requirement of near resonance. 

Perhaps more seriously, it is not obvious that there really is no orthogonality blocking in the ferritin. The intrinsic 
blocking will be very small - the only obvious other sources are either the dipolar interactions between the ferritin 
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molecules, "loose spins" on the surface of the ferritin, or strain fields in the sample. The effect of dipolar fields 
is probably small, but the other effects are difhcult to quantify - we simply call attention to the severe effect on 
the lineshape that will be wrought by even a fairly small amount of orthogonality blocking (see Fig. 16, for k — 2). 
Orthogonality blocking will also further reduce the spectral weight (Eq.(5.17)). On the other hand nuclear spin 
diffusion between the protons may be important. 

Thus there is apparently a disagreement between the theory and experiment for ferritin. Nevertheless it may not 
be insurmountable, provided the discrepancy concerning power absorption can be dealt with and provided 

one can rule out any significant orthogonality blocking, and we certainly do not think that one can yet rule out MQC 
as a possible explanation of the experimental results. Another question that needs further investigation is the role 
of the dipole interaction between the ferritin molecules - this may well lead to some kind of cooperative behaviour 
between the ferritin molecules. Elsewhere we have also suggested Q a search for a "spin echo" which would help to 
reveal the source of the experimental resonance at ~ 1 MHz. 

It should also be noted that very recently Garg has given an analysis Q| of the effects of the nuclear spins on 
coherence in ferritin molecules. Garg concludes (in agreement with our earlier papers that the experiments do 

not agree with the theory. However his reasoning and detai led r esults are different, and are based upon the reduction 
in the total spectral weight of the signal (the factor / in Eq.( 4.37 ), where / = (2/7riV)^/^; recall this factor comes solely 
from the assumption of zero nuclear spin polarization) . In effect the analysis of Garg ignores topological decoherence 
(afc = 0), ignores orthogonality blocking (k — 0) and it ignores degeneracy blocking {SuJk = 0). He therefore arrives 
at a set of sharp lines for the spectrum of x"(tt'). 

However, as we have just seen this is not realistic. Dropping ak is justified for ferritin (but not for most other 
systems - compare the values for the examples in section III.B), but we see that the spectrum in Fig. 13 is not a sharp 
Lorentzian line. Thus the analysis of Garg misses the crucial physical mechanism (of decoherence); his lineshape is 
one of completely coherent motion for the ferritin Neel vector (albeit with a reduced total signal amplitude) . We note 
moreover that the amplitude reduction factor found by Garg is just the "bare" value /, which is incorrect; the correct 
reduction factor is that given in Eq.(4.37), for pure degeneracy blocking (or (5.17), if orthogonality blocking is also 
present) . 

A further point, upon which we will elaborate in detail in another paper dealing with finite field effects, is that the 
field dependence of x"(a;) will be quite different from that of a 2-level system, even one coupled to a bath of oscillators. 
This is because degeneracy blocking causes x"(w) to vary on an energy scale determined by W{t — jeSHo),, which 
has nothing to do with Aq and is usually much greater. 

We hope that from our analysis of these 2 examples the reader will now have a fairly concrete idea of how our 
effective Hamiltonian is used in the analysis of a real central spin; in fact with the analytic results we have given, 
it should be possible for experimentalists to produce plots of x"(^) for given system, provided the relevant 
parameters are known. 



VII. CONCLUSIONS 

In this paper we have given an analysis of the behaviour of a central spin S, coupled in an essentially arbitrary way 
to a spin environment. That one can do this, and moreover recover closed analytic forms for the correlation functions 
Pfflf (t) and x"{^), is perhaps surprising. It was certainly not obvious to us when this work was initiated (cf. ref. |l|); 
although the 3 mechanisms involved in the decoherence of Pff{t) are already contained in the analysis of ref. [2], it 
has taken some time to find closed analytic forms for all cases. It is perhaps even more surprising when one examines 
the equation produced by the only other method we could think of using for this problem, viz., the Bethe ansatz (we 
are dealing with a set of TV coupled 1-dimensional problems here ]47| , with in general random couplings). However 
examination of the relevant Bethe ansatz equations reveals them to be apparently intractable. 

Our solution relies on our novel technique of introducing "operator instantons" in the Hilbert space of the envi- 
ronmental spins, and then separating out the relevant energy scales. This then allowed identification of 3 important 
limiting cases of the general problem, each of which involved a single specific physical mechanism. These mechanisms 
were " topological decoherence" (already completely solved in ref. ) , and " degeneracy blocking" and " orthogonality 
blocking" , whose essential physics was also described in ref. . The solution is then achieved by combining these 3 
mechanisms; this allows full coverage of the parameter space when there is decoherence without energy dissipation. 
In this paper we have used these methods to examine specifically the coherence properties of S in zero external field. 
This allows us to deal only with the zero polarization sector of the environment. In fact a fairly trivial extension of 
these methods can be used to deal with the finite field case, and to include all polarisation sectors. This then yields 
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a theory of the relaxation of S{t) (including tunneling of S{t)) in an applied field, or in its absence. This will be 
discussed in detail in another paper. 

As a result of these investigations we find that the spin environment acts very differently from the usual " oscillator 
bath" environments, and has a far more powerful suppressive effect on quantum coherence at low temperature. One 
can also analyze the effect on coherence, of spin environments, for other systems. For example, an analysis of possible 
macroscopic quantum coherence for SQUIDs yields quite enormous values of our parameters A, A', etc. (typically 
A > 10*^). This implies that the previous analyses of SQUIDs, with a view to observing such "MQC", are far too 
optimistic, and that in fact it is unlikely that it will be seen in any superconductor (although the prognosis is more 
hopeful if the superconductor lacks nuclear spins). This work will be reported elsewhere. 

These results lead us to suspect that in nature the overwhelmingly dominant mechanism of dccohercncc on the 
mesoscopic or macroscopic scales is likely to come from environmental spins, particularly nuclear spins. Whilst this 
may not change the fundamental attitude of most physicists towards the foundations of quantum mechanics, or the 
" measurement problem" , it is certainly an interesting and non-trivial discovery about the physics of decoherence in 
the real world. It is also of fundamental importance for efforts to either observe MQC, or to at least build devices 
which operate coherently on more than just the 1- or 2-particle level. As noted in the introduction, it very severely 
restricts the possible physical systems in which such coherence might be seen - the only really clear example that 
avoids the problem discussed here is superfluid ^He - again, the details of this example will be discussed elsewhere. 
Our prognosis for attempts to see MQC in magnetic systems is somewhat pessimistic, although if a system can be 
found with A ^ 1 and k <^ 1, we may be in business again. 

Finally, let us note again that once an applied field starts to noticeably bias the potential Ho{S), things change 
a great deal. Our problem then becomes effectively one of the effect of nuclear and other spins on the crossover to 
dissipative tunneling of S{t). This is a quite different problem from the one we have just analyzed, and we will return 
to it in another paper. 
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APPENDIX A: DERIVATIONS FOR THREE LIMITING CASES 



We describe here the derivations of the key formulae in section We begin with a few remarks on the free 
central spin, and the derivation of P-(jij{t) for it; it will then be seen how the expressions (4.4), (4.?;) involve obvious 
generalizations of the free spin algebra. Then we give an account of the lengthy derivation required to obtain the 
expression (4.22) for pure orthogonality blocking. 



1. Free central spin, plus topological and bias effects 

Once Ho{S) has been truncated, an expression for P^-f^{t) can be derived using standard instanton techniques |]l2| , ^ , 
employing the free amplitude = AqT^: exp{±i7rS'}. In the absence of the phase wS, one has a "2-level system" 
correlation fimction PrLsit) 

PTLsit) = ^[l + coa{2Aot)], (Al) 

which is derived in the instanton framework by summing over all paths involving an even number of flips, with an 
amplitude iA^dt to flip in a time dt, giving 

W)-2|l + E^^|. (A2) 

Adding the phase ttS clearly changes the flip amplitude in time dt to 2iAo exp{±i7rS'}, and adding the ± processes 
then gives 
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= ^ + cos(4Ao cos TTS)t] , (A3) 

which can, if one wishes, be derived [0j3^ by summing over all paths with an even number of flips, and summing over 
all combinations of clockwise and counterclockwise flips: 



(2.)! 



{2s-ny.nl 



^i*o(2s-2ri) 



(A4) 



where $o — ttS for the free spin. This is equation ( p^ ) of the text. 

The generalization to topological decoherence is now easy. We either observe that the interference factor (e""^ 



2 cos t:S in (|A3|) will now be generalized via 



N 



7r5 ^ ($ + ^ akfik ■ CTfc) ; (pure topological) , 
fc=i 



(A5) 



which gives Eq.(4.4) (including the trace over spin states); or we simply include the extra phase (A5) into (A4), to 
give Eq.(4^) of the text. 

Now consider the modification introduced by a bias e. For the 2-level system this is discussed very thoroughly in 
one way by Leggett et al. (ref. jl^. Appendix A). We use a different representation here. Consider first the amplitude 
A-ff-f^(t) for a return to an initial state in the biased case; this is clearly 



dhn ... dtie^(^(*-*-)-^(*--*--i)+--*i) , 
Jo 

±iedt 



since the action coming from a pause of length dt is e*""*. We now write this in Laplace transform form as 



A f \ 



(A6) 



(A7) 



p — le 



E 

n=0 



p2 ^_ g2 



2 \ " 



1 p2 



p — ie + E"^ 



(A8) 



which gives the standard answer for P^^{t^ e) 

dpidp2 



oo oo 

EE 



^ ' ^ ' n— rn— ^ 



(-*Ao)' 



A2 

1 %(1 -cos2M) 1 ^sin^ Et 

2^2 ^ ' E^ 



pl + e^ 



(A9) 



(AlO) 

It is then obvious that in the case of the free spin, for a bias that does not distinguish between the clockwise and 



counterclockwise flips, Pff{t, e) will look exactly like (AlO), except that Aq 2Ao cos nS. This representation avoids 
the difficulty in the usual representation which comes essentially from the square root E = [e^ + A^]^/^ in the cosine. 



2. Orthogonality blocking 



(1) 



We consider the situation described in section [V,B. For the spin ak, the "initial" and "final" local fields are 7^ 

and jj^'' ■ We will begin, to make the derivations easier, by assuming that the angle relating 7^^-' and 7^^-' is close to 
TT for all the at ■ This angle Pk is defined as 

cos2/?fc = -7f^7f^ (All) 
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We now choose axes such that the final state wave-function of Bk is given in terms of the initial state wave-function 



by 



(A12) 



JV 



{'?{}> = n I {'^n) = tj I {aD) 



(A13) 



fc=i 



In general the initial state of the spin bath will belong to the polarisation group AA'' (not necessarily AA'' = 0) , where 
the polarisation is defined along the S'l-direction. As explained in the text, when S rotates from to S2 ~ —Si, the 
energy conservation requires to change the polarisation from AN to — AN. In future discussion of degeneracy blocking 
effects we will have to consider more general tunneling processes, when the nuclear polarisation is changed from AA^ 
to AA^ — 2M and back. In what follows we calculate the correlation function PAN.M{t), which gives the dynamics of 
S when the spin bath transitions are restricted to be between the (P) = ASf and (V) — AN — 2M subspaces, which 
are supposed to be in resonance (for pure orthogonality blocking M — ZW). The statistical weight of states with 
AN ^ 7V^/^ is negligible, so we will assume in our calculation that AN, M <C N . 

The above restriction is enforced by the projection operator 



n 



ZW 



(A14) 



We can now, including this restriction, write down an expression for the amplitude for S to go from | Si{t = 0)) to 
I Si{t)) starting from some {V) = AN environmental state; 



A 



Kn=0 ^ ' i=l 



2tt' 



g iAiV(52„+e2>.-l+.--+Cl)g2iAf(e2„-l+«2„-3 + ---+?l)jl 



2« I {c?n> , 



(A15) 



where T2n is 



u 



From (A15) we can now write the correlation function PAN,M{t) as 

PAN,M{t)= (^AiV.Af (O^AiV,J\/(0) 

°° °° (lA t)2("+™) ^" ^™ 



EE 



(AI6) 



riff/ t/f ^-^^^^'"^-"^^"^^^^e-^^^— • (A17) 



(2n)!(2m)! 

n=Om=0 \ ' ^ ' i=lj=l 

The restrictions required to arrive at this formula are discussed in the text. 

We now return to the assumption that the Pk are small; specifically we will assume that the orthogonality exponent 
K, defined by 



N 



(A18) 



k=l 



can be approximated by the perturbative expansion 



1 ^ 



(AW) 



fc = l 



This assumption makes it much easier to calculate the average in (A17). Consider first the problem with only one 
environmental spin Ok, and calculate the average {T2rn^'^ri)k in this case; since T2„T2n is a product of operators acting 
separately on each dk, the average over all spins is also the product of single spin results. 
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Because of (A19) we need only to consider processes with 0, 1, or 2 flips of the environmental spin, i.e., we will 
make an expansion in powers of Pk for ak, and stop at Then it is clear that, if the initial state of ak is |tfc) 



(i-npi) \u)+^f3k iifc>E(-i)'+'e-2'n:,e 

1=1 

2n 2n-l 



(A20) 



where the first term arises from the sequence [ttt ■ ■ • TT]j the second from the sequence [ttt ■•• Tiii ••• ii]-, with a 
flip when j ~ I; and so on. In the same way we find 



2n 2n-l 



2m 2771—1 /_-|^ 



i'=i+i 1=1 



p'=p+i p=i 

2 m 2n 



p=i 1=1 



(A21) 



to order /3^. The sequence {lk\ {T2m)^'^2n life) ^il^ have a similar expression, but with reversed signs coming from 
the e^^^'^'k factors. 

We now observe that the state with polarisation AN consi sts of {N + AN)/2 spins up and {N — AN)/2 spins down. 
Consequently, for each ak, we add t or J, averages like ( A21 ), and then take the product 



Nr 



N, 



k=l k' = l 



Substituting (A21) into this expression we get 

{TLf2n) = e^'^T «-Er exp { - 

where the "effective action" K^{^{^i,£^j,AN) has two contributions K'^^^ = Ki + K2: 

I'-i p'-i 



(A22) 



(A23) 



AN ( , , ^ ~ 

K, = 2^c(l - — )| [n + m) + E(-l)' ' cos[2 E + E (^l)' '^"^P E ^'^ 

^ l'>l i — l p'>p j—p 

2m 2n 27i p—1 n 

-EE(-i^''-«[2Ee.-E^^]l 

p=i 1=1 i=i j=i ■' 



(A24) 



^2 = 2^^^^ E(-1)''"'°^P[-2*E^^] + E(-l)'''"''exp[2* E C,] 



N 



I'-i 



p'-i 



i'>i 

2m 2n 



i=l p'>p 
p-1 2n 



EE(-i)'^^-p[2^E^;-2zEe4 

p=l 1=1 j=l i=l 



(A25) 



We recall now that ZW < A^^/^ ^ N , which allows to neglect the contribution due to K2 and drop the correction 
AN/N to the coefficient k in Ki. Notice also that the phase factor in front of exp{— in (A22) cancels exactly 
the phase proportional to AA^ in the formular (A17) for PAN.Ai{t)- Thus, quite surprisingly, we find the correlation 
function to be independent of AN in this limit: 
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^mW-EE (tw2mV nn / £ / ^exp{2zM(6„-i+6„-3 + ...+Ci)-i^;i{{te,C;)}, (A26) 

n=Om=0 ^ ''^ 1=1 j=l-' •' 

We can render this expression more useful by changing variables; first we consider the whole sequence — 
(^1: • • ■ J S,2m ■ ■ • 5 ~C2m) together, and then define new angular variables 

2(n+m) 

Xa= E 2Ca + 7ra, (A27) 



so that now 



^mw^ 1:1: <||^ n /It 



exp|iAfE(-l)"+'Xa-2K[(n + m)+ E cos(xa - Xa')] [ 



(A28) 



Thus we have mapped our problem onto the partition function of a rather peculiar system of spins, interacting via 
infinite range forces, with interaction strength 2k. 

To deal with this partition function , we define "pseudo-spins" Sa and 5, via 

Sa = (cOSXa, sinXa) 
2(n+m) 

S= E ' (^29) 

a=l 

SO that 

Sa • Sa- = COs(Xa ~ Xa') 

Ecos(xa-Xo')=5^ , (A30) 



We can think of Sq as rotating in our fictitious angular space defined by the projection operator (A14). Now consider 
the term G{S) in (|A28[) defined by 

rfx«^.M(-i)-^x= I exp{-KEcos(xo-Xa')} 



(2(n+m) 
n / 
Q = l 

(n+m) 



This is easily calculated, viz. 



2(n+m) 



C( = l 



,^ „ / „2ii- 7 \ 2(n+m) 



27r 7 27r 

= ^/d-e-V4.4»+™)(,)^ (A32) 

where Jm{X) is the zeroth-order Bessel function. Using 

(2(n + m))! r,2(« + m)-l |^Ann^ 

(2m)!(2n)! "'^"+-" + 2 ' (A33) 
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to reorganize the sum over n and m in ( A28| ) and changing the integration variable z 2x^/k, we then find 



2 s 



2/ dxxe-"'p{y(i,AM(x)) . (A34) 







Am (a;) = AoJAf(2V^x) . (A35) 

Here we come to the crucial point in our derivation. Eq.(b.21b) gives us the final answer as a superposition of 
non-interacting correlation functions for effective tunneling rates Am (x) with the proper weighting 

Pni*) = dxxe~''\l + cos[2AoJMi'2VKx)t]) , (A36) 
Jo 

For M = this is the form quoted in Eq.(O^) of the text. 



It is worth noting that non-zero M enters this calculation as the overall phase factor starting from (A17) and crown 



it in (A32) while finally intergating over {xa} to produce the Besscl function of order M. This observation allows 
to generalise any calculation done for M = to finite M by simply replacing Jq Jm in the final answer - the 
prescription which we make use in other Appendices. 

APPENDIX B: DERIVATIONS FOR THE GENERIC CASE 

We outline here the derivations for section incorporating all three mechanisms. 

1. Orthogonality blocked topological decoherence 

As discussed in the text, a formal expression for this case can be easily written down, by including the induced 
topological phase ±[<i> + ^^^li'^k^k — i£,kVk) • o'fc] which accumulates during each instanton, and then averaging over 
the trajectories of the {at}- Here we evaluate the AI — contribution; The result for PM(t) then follows from obvious 
generalisation. We have 



Poit) -1^1^ (2n)!(2m)! ^ ' ' [[[[^ ^{TLim)TM) , (Bl) 

n=Om=0 ^ I {g,=±} i=\ 

where the transition matrices now include these phases, i.e., 

T2„(gi) = e'«-^i7/„pe^«--^^ . . . e^^^Utov ■ (B2) 
and f/top now generalizes the [/ involved in orthogonality blocking, viz.. 



AT 

Uto-o = exp \igi\(aunk - iikVk)\ ■ - j/3fcO-fc| . (B3) 



' top — 

fe=l 



To evaluate (Bl), we rewrite it as 

^oW - E E ^ t V(2^). 5: ^' n /f exp{-i.S(M,{U)}, (B4) 

n=Om=0 ^ ''^ {gi=±} p=l 

and calculate the K'^d ({5i}i{Cp}) np to order in a\, 0f,, £_l (see Appendix A. 2). 

Consider first the case where f3k and sue zero. Then, with the usual assumption that the individual ak are small, 



but not necessarily A = ^f^^^al, the "effective action" K^(^ simplifies to (compare Eq.(A28)) 
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2(n+m) 2(n+m) 

Kr/Ji{^p}) 3p9p' + (A - A') ^ cos(xp - XpOSpV , 

p' -.p p'^p 



(B5) 



which generahzes from orthogonahty blocking; the Xp £^re defined as in (A27), and 



k=l 
1 ^ 

A' = iE"'("^)' 



(B6) 



fc=i 

The term in A' comes from the expansion to 0{a^) (analogous to the expansion in (3f. in Appendix B), which 
produces the average 



(B7) 



(compare (A20); again, we assu me eq ual statistical weightings for t and [ spins). 
We now use the same trick of ( A29 ) and ( A30 ), to write 

2(n+m) \ ^ 

^ ^ \ X! ^P^P ^ Ysp9p' cos(xp - Xp') : 

P=l / P':P 



(B8) 



with the Sp = (cos^p, sin^p) as before; making the change of variables Xp — Xp + ""ffp/^, we get the generalization of 
(^H), viz., 



/If °° 
2 L „_n 



(2^)! 



{9l=±} 



2 



1+E 



[2iAoJo(2aVX^)t]2« 



s=0 
{9l=±} 



(B9) 



Writing the sum over {gi} in the form 

{9l=±} 



2(n + m) ^2 

91) _ 



2s 



(2s)! ij»(2s-2m)-A'(2a-2m)^ 



V — e 

^ (2s - m)!m! 



(BIO) 



we now use the identity in (4.8) to eliminate this term in favour of a phase integration over ip; defining _F\' (f^) = e 
as in the text, we get P^^(f) in the form 



where 

Ao(iy9,a::) = 2AoCosiy9Jo(2a;\/A — A') 
If we do the integration over Lp, this gives 



Po(t) = /te-"' |l+ £ FA'(m)e^2m$j^^[4^^^j^(2a;\/A^)]l , 

7n— — oo ) 



(Bll) 



(B12) 



(B13) 



35 



which is Eq.(5.3) in the text. 

The above calculation is generalized for non-zero M by noting that the only modification will be in a phase factor 
iM^p{—\Y^^Xp throughout the whole derivation, which will finally result in the replacement Jq{2x^/\ — A') 
JM{2xy/X — A') in Eq.(B12) for the effective tunneling rate. 



2. Orthogonality blocking plus degeneracy blocking 



We wish to evaluate the "biased orthogonality blocking" problem described in section We generalize the 

expression (A8) for a biased 2-level system (topological phase is not yet included here) to include orthogonality 
blocking via our projection operators (cf. Eq.(A15); this gives 



AM=o{p,e) 



(B14) 



with AM=o{t,£) given by the Laplace transform (A7) of ( B14 ). It then follows that 



oo oo 



dpidp2e 



{Pl+P2)t_ 



Pi -iep2- le 



^— X! X! Bnjn{pi,P2,e) 



(B15) 



n—0 m—O 



BnrniPl,P2,e) 



pl + e^ 



(»A^)^ 
pl + e^ 



1=1 j=i 



(B16) 



Surprisingly, we find that the problem of calculating the correlation function in a bias (after writing it as a series 
expansion in the tunneling rate) can be again solved in terms of non- interacting P^^{t,e), because the averages in 
( B16| ) are identical to those evaluated in Appendix A, Eq.(A32), to give 



2n 



2in 



27r 



(B17) 



This can be absorbed into defining an effective Ao(a;) — 2A$ Jo(2a;Y^), just as for pure orthogonality blocking, and 
the series become identical to Eq.( |A10 ) giving the result (5.10) in the text. The generalisation to M ^ just consists 
in the replacement Ao(x) — > Am (a;). 



3. The generic case 



parallel 



, <&, and ^fc) except uj-f. 



perp 



included from the effective Hamil- 



Here we calculate P^^(i) with all terms (ok, 
tonian. Thus we return to the expression (|B3| ) and (B4) and calculate the "effective action" K^{yl with nonzero ^fe. 
Now , evaluating the contribution from individual spins from the expansion to 0(a^., d), we find for the average (cf. 

(H) 



X {iakfik + [-ly ^kVk) ■ CTfc I CTfc) , 

the following expression 

(lakul + i-ir^kv'k) {^ockul + {-ly'ikvl) 

, -2ty'-u, 



(B18) 



(B19) 



Thus, in full analogy to the derivation done in Appendix B (Eq.(A21)), we find the contribution of an individual spin 
to the action to be 
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2(n+m) , 



- lakUnk ■ Vk - nlvm-^r + (-1)"')] } , 



(B20) 



Summing up the contributions from spins we generate the "effective action", which now takes the form 

2(n+m) , 

p,p'=i 

+ cos(xp - xp') [(A - A') - i-ir+p'iv - v') - - m-^r + i-ir')] } , (B2i) 

where the constants are defined by: 

N N 

A=^E"^ A' = iE"^("^)'; (B22) 



2^"' 2 

k=l k=l 



N N 



2 ' ' 2 

fc=l k=l 



N , N 



7 = ^ E ^'k^krik • 4 ; 7' = ^ E '^kikulvl ; (B24) 

k=l k=l 

As before we change variables according to Xp = Xp + tt when gp = — 1, to introduce the odd and even spin fields 

2(ri+m)-l 2(n+m) 
p=odd p=even 

Now we notice that the effective action is a quadratic form in (5o,5e): 

ifn{f = A'(E5p)'-^'(E(-i)%)'-2n'(E5p)(E(-i)%) 

p p p p 

+ (A - A')(5; + 50)2 - (r^ - - 5;)^ - 2z(7 - 7')(52 - , (B26) 

2(n+m)-l ^ N 
p=odd -^pi 



P P P P 

^2 „'^/'c c a2 r,„-/'_, -/\/'c2 c2^ 

The rest is simple now. First we employ spectral representations for the (5-functions 6{So — X^p^odT^ ^ ^p) 



S{Se — X]p=euen ^p) integrate over the angle variables Xp (compare Eq.( A32 )) to give (s = 71 + 111) 



p=l"' ^ 

X exp j - a{Se + So)^ + b{Se - Sof + 2ic{Sl -Sl)+ixi-Se + ix2 ■ So\ ■ (B27) 



with obvious definition 

a = A-A'; b = t] - r]' ; = 7-7'. (B28) 
The Gaussian integration over dSo, dSe, and then the integration over the angle between the xi and X2 gives for the 



expression (B27) the formula 
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1 / (a + h)xiX2 \ f (a — 6 + 2ic)x\ + (a — 6 — 2ic)x\ 



8(a6-c2) " V 8(a6-c2) ; ^[ 16(a6-c2) 
Finally we note that the sum over {gp\ can be written as 

oo 

E ^ E E mmi-Y.9p)5{2m2-Y.{-lYgp) 

{gp=±) mi,m2 = -oo {gp=±} p p 



(27r)2 

mi,m2 = — oo {gp=±} 

E e^2™,^,+.2™2^2 [2eos(^i + 2cos(^i - ^2)] ■ (B30) 

^ 7ni ,m2 — — 00 



Combining now (^), (^ ), and (HI) together we observe that the sum over n and m is nothing but the series for 



the coherent dynamics expansion in powers of the renormalized tunneling splitting 

Al{xi,X2,ipi,ip2) =4A^cos(i^i +(^2)cos(v3i - ip2)Jo{xi)Jo{x2) , (B31) 
which is then integrated over the variables xi,X2,y^i, f2 and summed over toi, m2 with the weight given by 

2 — g2i[mi(<I> — ipi)— m2¥'2+4mim27']g4(i)'m| — A'mj) 

xix-? (fl + b)xiX2 \ f {a — b + 2ic)x\ + (a — 6 — 2ic)x\ 

One then finds 



''W^)'"\ i{ab-c') j""! 16(at-c') 



•' •' mi=-oo m2=-oo 

X Z(.^i,(/92,a;i,X2,mi,TO2)P]:',|-'[t, Ao(a;i,X2,(/?i,(i52)] , (B33) 



and (B33) has an obvious generalisation to include the bias integration / de. 

FIGURE CAPTIONS 

Figure 1 Schematic depiction of the Central Spin Model 

Figure 2 The range of coupling energies between a typical central spin and its spin environment. The hyperfine 
contact interaction is shown for a variety of nuclear spins in various hosts. 

Figure 3 A typical trajectory for S'(r), contributing to the path integral. 

Figure 4 The function OLoiy^o) for the instanton of the toy model Hamiltonian. This form is typical of virtually 
any instanton describing the tunneling of our central spin. 

Figure 5 Behaviour of the central spin in the topological decoherence limit. In (a) is shown P-^-^(i) for the 
intermediate coupling case; and (b) shows the spectral function x"('^) deriving from this. These forms are the 
universal forms arising for topological decoherence (see text) . 

Figure 6 P-f[f{t) for an 5* = 1/2-integer spin coupled to 40 environmental spins with random ak between and 
0.05. 

Figure 7 Pffit) for an integer S, and coupling to 20 environmental spins with random ak between and 1. 
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Figure 8 Definition of the orthogonality angle in terms of the unit vectors 7^^-' and ■ This also defines the 
unit vectors Ik and rfik appearing in the effective Hamiltonian (3.10). 

Figure 9 The effect of pure orthogonality blocking when n is small (k < 0(1)). 

Figure 10 The ef fect o f pure orthogonality blocking when k ^ 1. As well as being pushed to lower frequencies in 
accordance with Eg. ( 4. 27 ), \"{uj) develops a complicated peak structure. 

Figure 11 P^-^{t), in the orthogonality blocking limit, for n — 1000. Notice the small resurgence of oscillation for 
Aot ~ 300, coming from constructive interference between the environment-induced peaks in x" {lo). 

Figure 12 The effect of a finite but small spread in the couplings {wfc} on the energy levels of the system. In (a) 
the set of levels for the total "central spin + environment" is shown as we add environmental spins one by one. In (b) 
the resulting distribution W{e) is shown. The parameter ^, = N^^'^Sujk/wo- The separation between peaks for /i ^ 1 
is iOo- The tunneling splitting is far smaller, indeed invisible on this scale for typical values of Ug- Usually > 1, 
giving the smooth Gaussian; if /i <C 1 a peak structure appears in W{e). 

Figure 13 The spectral weight x"{^) for the case of pure degeneracy blocking. 



Figure 14 A plot of the result for projected topological decoherence, in Eq.(5^), for small values of (A — A'); 
notice that even if (A — A') is only 0.1, the result is very noticeably different from pure topological decoherence (for 
which A — A' = 0), shown in Fig. 5b. 

Figure 15 Graphs of x"{^) foi' projected topological decoherence, but now for A — A' 1. 



Figure 16 The spectral function x"(w) for degeneracy blocked orthogonality blocking (eq.( 5.1(: )). Here we see 
the results for small n . 

Figure 17 x"{^) for degeneracy blocked orthogonality blocking, for 1. 

Figure 18 x"{^) for the case where all 3 decoherence mechanisms enter, and we assume that A, A' > 1 (so there 

is topological decoherence). Here we show x"(<^) for small values of A — A'. 

Figure 19 x"{^) when all 3 mechanisms operate, and A, A' > 1; here we show results for large (A — A'). 

Figure 20 x"{^) for ThFe^ grains, calculated using the numbers in the text. 
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TABLE I. SELECTED HYPERFINE COUPLINGS 



Nucleus 


Host 


State 


Wfc (MHz) 


Abundance (%) 




CuCh ■ H2O 


(AFM) 


3.5 


99.99 


pl9 


MnF2 


(AFM) 


160 


100 




MnF2 


(AFM) 


680 


100 




MnB 


(FM) 


217.7 






Fe 


(FM) 


46.65 


2.19 




a — Fe20z 


(FM) 


71.5 






FesOi 


(Fcrritc) 


63.55 


75 




YiFesOi2 


(YIG-FM) 


64.9/76.05 


75 




FeCo5; 


(FM) 


44.8 


75 




FeCo5 


(FM) 


289.2 


100 




Co 


(FM) 


228 


75 




Ni 


(FM) 


28.35 


1.19 




Tb 


(FM) 


3776, 3108, 2439 


100 




Dy 


(FM) 


1603, 1219, 830, 445 


18.88 








1984, 1574, 1163 


24.97 








755, 347 




TABLE 11. 


SOME TYPICAL NMR LINEWIDTHS IN ANTIFERROMAGNETS COMING FROM THE SUHL - NAKA- 


MURA INTERACTION. For Mn^^ a 1 Oe linewidth roughly corresponds to a Jwfc ~ 1 kHz. 




Nucleus 


Host 






Linewidth (Oe) 




MnF2 




(AFM) 


14 








(AFM) 


1000 




C0F2 




(AFM) 


268 (1/2 < — » -1/2) 








)7 


258 (±3/2 ^ ±1/2) 




55 




)7 


222 (5/2 < — > 3/2) 


): 


7) 




)5 


165 (7/2 ^ 5/2) 
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